MUTUALLY CATALYTIC BRANCHING PROCESSES AND VOTER PROCESSES 

WITH STRENGTH OF OPINION 



LEIF DORING AND LEONID MYTNIK 



Abstract. Since the seminal work of Dawson and Perkins, mutually catalytic versions of super- 
processes have been studied frequently. In this article we combine two approaches extending their 
ideas: the approach of adding correlations to the driving noise of the system is combined with the 
approach of obtaining new processes by letting the branching rate tend to infinity. The processes are 
considered on a countable site space. 

We introduce infinite rate symbiotic branching processes which surprisingly can be interpreted as 
generalized voter processes with additional strength of opinions. Since many of the arguments go 
along the lines of known proofs this article is written in the style of a review article. 



Going back to the seminal work of Watanabe |W68| and Dawson |D78| . the subject of measure- 
valued diffusion processes arising as scaling limits of branching particle systems has attracted the 
interest of many probabilists. Many tools had to be developed to study the fascinating properties of 
the Dawson/ Watanabe process (also called superprocess or super-Brownian motion) and its relatives. 
Characterizations and constructions of the process via a Laplace transform duality to non-linear 
parabolic partial differential equations, infinitesimal generator and corresponding martingale problem, 
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the pathwise lookdown construction of Donelly/Kurtz [DKlj . |DK2j or Le Gall's Brownian snake 
construction based on the Ray-Knight theorems (see the overview |LG99j ') led to many deep results. 
Much of the analysis is based on the branching property, i.e. the sum of two independent super- 
Brownian motions X} , X^ is equal in distribution to a single super-Brownian motion started at Xq + 
XI 

In the early 90s further directions became popular. Super-Brownian motion was found to be a universal 
scaling limit not only of branching systems but also of interacting particle systems such as voter process 
and its modifications (see for instance |CDPOO| . |CP05) ). Furthermore, instead of considering plainly 
super-Brownian motion, interactions were introduced. Tools such as Dawson's generalized Girsanov 
theorem |D78| have been successfully applied in various contexts. Here, we will be mostly interested in 
variants of catalytic super-Brownian motion, i.e. super-Brownian motion with underlying branching 
mechanism depending on a catalytic random environment. As long as the environment is fixed, a 
good deal of the analysis can still be performed with techniques developed for the super-Brownian 
motion. More delicately, taking into account connections to stochastic heat equations, Dawson/Perkins 
introduced a mutually catalytic superprocess (see |DP98| V Their mutually catalytic branching model 
on the continuous site space M consists of two super-Brownian motions each being the catalyst of 
branching for the other. The model was described via stochastic heat equations. They considered 

\§iut(x) = \l\ut{x) + ^-tUt{x)vt{x)Wl{dt,dx), t > 0,a; G M, 
I wM^) = \^vt{x) + ^-iUt{x)vt{x)W^{dt, dx), i > 0, X € M, 

driven by two independent white noises W^, on M+ xM. Here, A denotes the one dimensional Lapla- 
cian. The mutually catalytic interaction of two super-Brownian motions has one particular drawback: 
the branching property is destroyed so that many of the previously known tools collapse. Fortunately, 
some ideas borrowed from the study of interacting particle systems and interacting diffusion models 
could be applied successfully due to the symmetric nature of the model. In particular, a self-duality 
that extends the linear system duality known for interacting particle systems could be established 
and utilized to prove uniqueness and longterm properties. Besides the above continuous model, the 
mutually catalytic model on the lattice was constructed and studied by Dawson and Perkins as well. 

This article, which is focused on spatial branching models on discrete space, is motivated by two re- 
cent developments. First, in the series of papers [ KMlOa] . |KM10bj . |KM11-|-| the effect of sending the 
branching rate 7 to infinity was studied in the discrete space mutually catalytic branching model. The 
resulting infinite rate mutually catalytic branching model is one of the rare tractable spatial models 
with finite 2 — e moments but infinite 2nd moment forcing the system to have critical scaling behavior. 
Secondly, Etheridge/Fleischmann |EF04) introduced the following generalization of the Dawson-Perkins 
model. They considered the mutually catalytic branching model with correlated driving noises which, 
on the level of a branching system approximation, corresponds to a two type system of branching 
particles with correlated branching mechanism. They called their model symbiotic branching model 
in contrast to the mutually catalytic branching model of Dawson/Perkins that appears as a special 
case for zero correlations. We will use equally the name symbiotic branching and mutually catalytic 
branching with correlations. Correlating the branching mechanism might seem artificial on first view. 
On second view one observes that the extremal correlations lead to well-known models from the the- 
ory of interacting diffusion models: the stepping stone model with applications in theoretical biology 
and a parabolic Anderson model with applications in statistical physics. As those models have very 
different path behavior one could expect phase-transitions occurring when changing the correlations. 
On the level of moments those phase transitions have been revealed recently in |BDEllj : there is a 
precise transition for 2nd moments when the correlation parameter changes from negative to positive. 

The main result of this article, formulated here in a slightly simplified version, is the following theorem 
which should be viewed as the natural combination of the two aforementioned developments. In 
particular, the theorem below extends results from |KM11-|-| to the case of "correlated (symbiotic) 
branching" . 
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Theorem 0.1. Suppose g € (—1,1) is a parameter and [u] ^v^)t>o is the unique non-negative weak 
solution to the symbiotic branching model on the lattice defined by 



(0.1) 



'dut{k) = Aut{k) dt + y^-fvt{k)ut{k) dBl{k), fc e Z'^, 



\dvtik) = Avt{k) dt + ^-ivt{k)ut{k) dBf{k), keZ'^. 
Here, A denotes the discrete Laplacian on 1/ 

\i-k\ = l 

and the driving Gaussian process {B^ (k) , (k)} f-^^d has correlation structure 

E[Blik)Blij)]=Soik-j)t, 
(0.2) E[BUk)Bf{j)]^Soik~j)t, 

E[Bl{k)B^{j)]^gSo{k-j)t. 

Additionally, assume that the non-negative initial conditions (mq,Wq) = {Uq^Vq) do not depend on 7, 
satisfy a minor growth condition (for the precise definitions see (S.S^ and Section \2.1.1 ) and also 

{Uoik),Voik)) eE:^ {(yi, 0) : yi > O} U {(0, 2/2) : 2/2 > 0} C 

for all keZ'^. 

Then [u'^ ,v'^) converges, as 7 tends to infinity, weakly in the Meyer-Zheng "pseudo-path" topology 
(introduced in [MZ84] ). to a limiting ROLL process {Ut,Vt)t>a taking values in E which is the unique 
weak solution to the system of Poissonian integral equations 

Utik)^Uo{k)+ f AUs{k)ds 







■ {y2V,-{k) + (2/1 - l)[/,_ (fc)) {N - N'){{k}, diy^,y2),dr, ds), 



Jo 



(0.3) 

Vtik) = Vnik)+ I AVs{k)ds 





ft fhik) j- 

{y2Us-{k) + (j/i - l)Vs- (fc)) (AA - AA')({fc}, d{y,,y2),dr, ds), 

lO Jo JE 

for t > and k G Z"*. Here, M is a Poisson point measure onl/ x E y. (0, 00) x (0, 00) with intensity 
measure 

Af'{{k},d{yi,y2),dr,ds) ^ iy<^(d{yi,y2))drds, Vfc G 

du : 2/2 = 0, 

^^(%i,y2)) ■ ^' 
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V^t-(fc) >0, 
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Remark 0.2. Theorem \0. l\ will be proved in Section^for more general countable state -spac e S instea d 



ojTL and Q-matrix A instead of A. The proof of the theorem follows from Theorems 3.4 



and 



3.15 



The parameter g only occurs in the measure so that it does not surprise that proofs go along the 
lines of IKMll+1 replacing in their Poissonian equations p(0) = 2 by some p{g) > 1. The striking fact 
of the generalization to p 7^ is that it allows to understand {U, V) as a family of generalized voter 
processes with the standard voter process appearing for g = —1. 

The generalized voter process interpretation goes as follows: Suppose at each site k E Z'^ lives a 
voter with one of two possible opinions. Their opinions additionally have a non-negative strength. 
Mathematically speaking, the type of opinion is determined by the non-zero coordinate of the opinion- 
vector (recall the definition of E) and the strength is determined by the absolute value, i.e. 

• (u, 0) E E codes opinion 1 of strength u, 

• (0, v) & E codes opinion 2 of strength v. 

Formulated like this, the standard voter process only takes values (1,0) and (0, 1) since all opinions 
do have a fixed strength, say 1. If u (resp. v) is large, we say the opinion is strong, otherwise weak. 
Voters change dynamically their opinions and their strength according to the next two possibilities: 

• Change of opinion strength only: Suppose M has an atom at (A:, (yi, 0), s, r). Then, by 
definition of the two integrands, the Poissonian integrals produce two-dimensional jumps of 
the form 

so that, added to the current state of the system, the state of the system at site k changes 
according to 

'Us-{k)\ (Us-{ky 



If before the jump the voter had opinion 1 of strength m, the change is (u,0) 1— > (yiit, 0) and 
(0, v) I— (0, yiv) if the voter had opinion 2 before. Hence, if (yi, 0) is chosen by the basic jump 
measure ^ only the strength of the opinion changes but not the type. 

Change of opinion and its strength: Suppose N has an atom at (fc, (0, t/2), s, r). Then, by 
definition of the integrands, the Poissonian integrals produce jumps of the form 

so that, added to the current state of the system, the state of the system at site k changes 
according to 

Vs-{k)) '^^\Us-{k) 

If before the jump the voter had opinion 1 of strength m, the change is (u,0) i-> (0, 1/2") and 
(0, v) i~7> {y2V, 0) if the voter had opinion 2 before. Hence, if (0, 2/2) is chosen by the basic jump 
measure ^ the voter changes strength and type of opinion. 



Remark 0.3. We show in Section \ 3.5\ that Theorem \ 0.1\ extends naturally to p = — 1 when is 
replaced by i5(o,i)- If additionally (C/o,Vo) G {(0, 1), (1, 0)}^ , then solutions to (0.3) give standard 
voter processes. Note that in this case only the second type of changes occurs since J\l only has atoms 
at (k, (0,1), s,r) . Hence, the strength of the opinion does not change. In particular, we only see opinion 
changes from (1, 0) to (0, 1) and vice versa. 

Finally, we should also give an interpretation to the rates It{k): due to the definition of It{k) and A, 
the rate of change for the voter at site k is high if the strength of the opinions of his neighbors of 
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different opinion is high compared to his opinion. In particular, voters with weak conviction tend to 
change quicker their opinions than voters with strong conviction. 



The resuh of Theorem |0 . 1 1 might look frightening to the reader not familiar with interacting diffusion 
processes and/or jump diffusions. However, once the connection to the results of |KM11+| and jBDEllj 
is understood, the proofs of the theorem go along the lines of IKMll+1 . Therefore, we decided to write 
this article in the form of a review article explaining in depth the background. We do not give many 



detailed proofs but instead give more detailed calculations to explain the origins of (0.3). In the 
following we explain carefully 

• the background of catalytic branching processes, 



definitions, existence, uniqueness and tools for (0.1) 



what is known on the longtime behavior of (0.1) to motivate the choice of in the theorem 
via planar Brownian motions exiting a cone. 



• more details on (0.3) and (alternative) constructions of {U,V), 

• concepts and definitions for jump diffusions. 

The background and connections to well-known stochastic processes from the literature will be ex- 
plained exhaustively in Section [T] Two different routes from known models to mutually catalytic 
branching models are disclosed: the original motivation of Dawson/Perkins originating from catalytic 
super-Brownian motion and symbiotic branching as unifying model for some interacting diffusions. As 
a final motivation the connection of stepping stone processes and voter processes is recalled. Section [2] 
is devoted to an overview of precise definitions, existence and uniqueness results and longtime prop- 
erties for finite rate symbiotic branching processes. In particular, the second moment transitions are 
discussed in detail. Proofs are cooked down to the main ingredients. Finally, in Section [3] the infinite 
rate symbiotic branching processes are introduced and reinterpreted as generalized voter processes in 
the very end. Additionally, a brief summary of jump diffusions is included to the appendix. 



1. Background and Motivation 

1.1. Prom Superprocesses to Mutually Catalytic Branching. Being a major subject of prob- 
ability theory, measure-valued diffusions, or superprocesses, such as super-Brownian motion and the 
Fleming- Viot process have been well studied during the last three decades. Important properties of 
superprocesses have been proved and connections to other areas of mathematics such as partial differ- 
ential equations have been established. For a detailed exposition of the subject the reader is referred 
to [D91], [Fog and EOO]. 

Here we introduce briefly super-random walks - the spatially discrete analogues of super-Brownian 
motion. Studying these processes gave a strong motivation to investigate spatial branching processes 
with interactions, and in particular, mutually catalytic branching processes on discrete space - the main 
theme of this article. To introduce super-random walks, we start with the following approximating 
particle system. Assume that an initial configuration of a large number (of order N) of particles 
distributed over lA is given. The particles move as independent simple random walk in 1/ and each 
particle independently of the others dies after an exponential time of rate , with 7 > 0, and at the 
place of death it leaves a random number of offspring particles, drawn from a fixed integer valued law 
jjL. The particles of the updated population continue their motion and reproduction according to the 
same rules. This process is usually referred to as a branching random walk with the branching law ^ 
and we will assume in the sequel that /i has expectation 1 (this means criticality) and finite variance 
cr^ > 0. The process x[^^ is then defined to be the finite atomic measure which loosely speaking gives 
measure of mass 1 to each particle alive at time t. To be more precise. 



i 
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where Xi is a position of the i-th particle ahve at time t. Assume that, as N tends to infinity, j^x'^^'^ 
converges weakly in the space of finite measures on to a measure Xq. Then one can show that the 
measure- valued process {-^Xj^-'}t>o converges weakly to a limiting measure- valued process {Xt}t>o 
which is called super-random walk and is uniquely characterized via the following martingale problem: 
for bounded test-functions : Z** — > M 

Mt{(j)):= f f lA(b{x)Xsidx)ds- f (t,{x)Xoidx) 
Jo JZ'' ^ Jta 

is a square-integrable martingale with quadratic variation process (7^7 /J /^d (f>'^{x)Xs{dx)ds. Here, A 
denotes the discrete Laplace operator as defined in Theorem |0.1| An interesting observation is the 
following invariance property: irrespectively of /i, the finite variance assumption for the branching 
mechanism ^ leads to a universal limit depending only on the variance and the parameter 7 which 
is also called the branching rate. In what follows, we assume cr^ = 1. It is worth mentioning that 
if we ignore the spatial motion and count just the total number of particles, the scaling procedure 
is nothing else but the scaling of critical and finite variance Galton- Watson processes which leads 
towards classical Feller's branching diffusion 



where Zt = Xt(Z'^), t > 0. 

Note that super-random walks can be characterized as solutions to stochastic differential equations. 
Abbreviating Ut{k) = Xt{{k}), the super-random walk is a weak solution to following system of 
stochastic differential equations (which is, in fact, a discrete version of a stochastic heat equation) 



(1.1) dut{k) = Aut{k) dt + y^-yutik) dBt{k), keZ'^, 

where {B{k)}f^^^d is a collection of independent Brownian motions. Next, we proceed to a more recent 
development: measure-valued processes with interactions. One way to introduce interaction into the 
model is to replace the constant branching rate 7 in the particle approximation by a random, adapted 
and space-time varying branching rate ^{t,k,uj), also called the catalyst. Some particular choices 
of branching environments 7 and related models over continuous space have been discussed in the 
literature (see for instance |DF94j . |DF95| . [D96j). For exampl e, on e can consider a super-random 



walk on 1/ in a super-random walk environment. Building upon ( 1.1 ), this model can be described as 



a solution to the following system of stochastic differential equations: 
(1.2) 

driven by independent families of independent Brownian motions. A solution (ut)f>o is called super- 



' dut{k) = Awt(fc) dt + ^vt{k)ut{k) dB]{k), k e 1"^, 
dvt{k) = Avt{k) dt + ^/-ivt(k) dB^{k), keZ'^, 



random walk in the catalytic super-random walk environment 7(-,fc,a;) = v.{k)(u)). Note that (1.2 1 



describes the so-called one-way interaction model: the v-population catalyzes the w-population. Then 



the natural extension of (1.2) to two-way interaction is the following mutually catalytic model. 



Definition 1.1. In the following, weak solutions {ut,Vt)t>Oj on a stochastic basis {^l^T, {J-t)t>Oi^), to 
the infinite system of stochastic differential equations (0.1^ driven by independent Brownian motions 
will be called mutually catalytic branching processes with initial conditions uo,vq and branching rate 
7 > 0. To abbreviate, solutions will be denoted by MCB^. In the sequel MCB^ will also denote a 
mutually catalytic branching process defined on a more general state space S instead of 1^ and with 
Q-matrix A instead of A. 

It is easy to see that the branching property fails for MCB^. Hence, many of the classical tools 
developed for superprocesses also fail. Nonetheless, the simple symmetric choice of the interaction 
between u and v makes this mutually catalytic system tractable. 
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Convention 1.2. In order to stress the underlying branching processes, the two components will be 
called types. 

As an example for the convention, if Ut{k) = for all fc e Z"^ wc will say that the first type died out. 

1.2. Prom Interacting Diffusions to Symbiotic Branching. Interestingly, the study of mutually 
catalytic branching processes can also be motivated by the study of interacting diffusion processes. 
Given a family of independent Brownian motions {St(A:)}^,g2d and some function / to be specified 
below, discrete-space parabolic stochastic partial differential equations 



(1.3) 



dwtik) = Awt{k) dt + ^-ff{wt{k))dBt{k), 
wo{k) > 0, ke 1'^, 



have been studied extensively in the literature. Some prominent examples will be briefly discussed in 
the sequel. 



Example 1.3. For J{x) — x solutions of are super-random walks. 

This example has already been dealt with in detail in the previous subsection. 



Example 1.4. For f{x) = x{l — x), Equation (1.3) is called stepping stone model. 



In fact, the stepping stone model is the spatial generalization of the one-dimensional Wright-Fisher 
diffusion 

(1.4) dXt = ^-iXt{l-Xt)dBt 

that arises as a scaling limit of the Moran model in population genetics similarly as the Feller diffusion 
arises as a scaling limit of critical Gallon- Watson processes. In contrast to the Gallon- Watson model, 
the Moran model is not used to model the total number of individuals but instead counts the proportion 
of one allele in a diploid population for a fixed number of individuals. In particular, this interpretation 



corresponds to the solution of (1.4) taking values in [0,1] with absorption at 1 or interpreted as 



fixation of genetic types. For an introduction to the questions of mathematical population genetics 



we refer to the lecture notes [E09| . The stepping stone model of Example 1.4 can be seen as an island 
version of the Wright-Fisher diffusion, i.e. additionally to the change of alleles, individuals live on 
islands which they change according to a nearest neighbor random walk. 

Changing the scope once more, we have a look at statistical physics. Given a random field rit{k), 
possibly time-inhomogeneous, the discrete heat equation with random potential 77 



(1.5) 



f ut(fc) = Autik) + i]t{k)ut{k), 
uo{k) > 0, fee Z'', 



has attracted a lot of interest. It is usually referred to as a parabolic Anderson model. Again, there is a 
connection to a branching particle system: started at localized initial condition uq — l{o}j Ut{k) is the 
expected number of particles in the system where one particle starts at and branches binary according 
to the breeding potential 77. In particular in the case of time-independent iid random potential a 
detailed analysis of the behavior of solutions is possible; we refer to the overview article |GK05| . If rj 



is the white noise case, then (1.51 is a particular case of (1.3) leading us to the next example 



Example 1.5. For f(x) = x^ , Equation {1.3) describes the parabolic Anderson model with Brownian 
potential (white noise potential). 

A detailed analytic study of the longtime behavior for this model can be found in the monograph 
|CM94| ■ For the probabilistic approach based on an explicit Feynman-Kac representation we refer to 
|GdH07] and references therein. 

Finally, the simplest example should be mentioned. Already in this case, a non-trivial interplay of 
noise and drift can be observed (see |CKOO] ). 
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Example 1.6. Choosing f{x) = I/7, Equation (1.3) describes interacting Brownian motions 



Now, as we have discussed examples that are of very different nature in terms of their origins and also 
of their properties, we should explain the connections to mutually catalytic models. Here is a prelimi- 
nary definition for the two types interacting diffusion model introduced by Etheridge/Fleischmann in 
|EF04] ■ A more precise and more general definition is given in Section 2. 

Definition 1.7. In the following, weak solutions {ut,Vt)t>o on a stochastic basis (fi, J^, (J^t)t>o, F) 



to the infinite system of stochastic differential equations defined in (0.1) driven by Brownian motions 



with correlation structure (0.2) are called symbiotic branching processes with initial conditions uqi^^Oi 



branching rate 7 > and correlation g £ [—1, 1]. 



To abbreviate, the system of equations (0.1) and their solutions will be denoted by SBM^{g) or just 



SBM^. 

The name symbiotic branching model was used in jEF04| in order to stress the biological interpreta- 
tion of the mutually catalytic behavior; the solution processes Ut and Vt might be considered as the 
distribution in space of two types. 

Convention 1.8. For later use let us capture the correlation structure used for symbiotic branching 
in a name. We will say that two Brownian motions satisfying E[i?j_B^] — gt, t > 0, are g-correlated. 

Having introduced the basic equations of this article, their relevance is emphasized by the following 
observation due to |EF04j . For correlation g = 0, solutions of the symbiotic branching model are 
solutions of the mutually catalytic branching model MCB-^. 

The case g = — \ with the additional assumption mq -l-wo = 1 corresponds to the stepping stone model. 
To see this, observe that in the perfectly negatively correlated case B^{i) — —B'^{i) which implies 
that the sum u-'r v solves a discrete heat equation and with the further assumption uo + vq = 1 stays 
constant for all time. Hence, for alH > 0, Wt = 1 — Vt which shows that m is a solution of the stepping 
stone model with initial condition uq and w is a solution with initial condition tig- 
Finally, suppose w is a solution of the parabolic Anderson model, then, for = 1, the pair {u,v) := 
(w, w) is a solution of the symbiotic branching model with initial conditions uq = vq = wq as now 
B^{i) = B\{). 

1.3. Infinite Rate Symbiotic Branching Processes and Voter Processes I. To motivate the 
procedure of sending 7 to infinity in Theorem |0.1| and to highlight for a first time why the generalized 
voter processes appear as limits, let us briefly discuss the voter process and its connection to the 
stepping stone model. For extensive information about interacting particle systems we refer to the 
monograph of Liggett |L05) . 

A way of defining interacting particle systems is a description via infinitesimal generators. Here, we 
assume that the voters live on Z"* and communicate only with their nearest neighbors. To define the 
dynamics via a generator, the state-space 

S = {0,lf' 

is fixed. The generator acts via 

(1-6) M^)= E cik,rj){firj('^)-f{ij)), 

on test-functions / : S ^ M only depending on finitely many coordinates and 77*^*^^ is defined to be 
the configuration in which the opinion is flipped only at site k and the rate of change at site k is 
proportional (normalized to total rate 1) to the number of neighbors with different opinion: 



2d ^ -{'7«#'/(fc)}- 

\i~k\=\ 
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Interestingly, the analysis of the longtime behavior of a voter process is drastically simplified by a 
pathwise graphical construction (see ^D93j . page 129): for each voter a vertical line is drawn downwards 
and each line carries a Poisson process firing tacks on that line. At each tack, a horizontal line is drawn 
randomly to a neighbor. With an initial configuration Uq G S, the construction goes as follows: for each 
site k with C/o(fc) — 1 water is filled into the vertical line and disperses downwards. Whenever there is 
an arrow pointing away from the line (this corresponds to persuading a neighbor) the water goes on 
downwards and, additionally, flows through the arrow to disperse downwards in the neighbor's line. 
When an arrow points from a neighbor's line towards the voter's line the water stops (this corresponds 
to be persuaded by a neighbor). At time t > 0, the configuration Ut is defined as follows: all sites 
filled by water carry a 1 and al 1 others a 0. It is heuristically clear that this construction yields a 



Markov process with generator (1.6) but interestingly it simultaneously gives a useful dual relation: 
reversing time and using the same arrows in the opposite direction, the resulting process is a system 
of instantaneously coalescing random walks. A simple consequence of this construction is a moment 
formula for the voter process: 

(1-7) E^'^[Ut{ki)---Ut{k^)] =E[[]c/o(ej)|eo'-fcir-- = 

i 

where ..-jC™ independent simple random walks that coalesce instantaneously when colliding. 
The product runs over all non-coalesced random walks at time t. 

Now, let us return to the stepping stone model 

(1.8) dwtik) = Awtik) dt + v/7Wt(fc)(l - wtik)) dBtik) 

that was already identified to the symbiotic branching process s g = —I. Unfortunately, there is no 
direct graphical construction for the stepping stone model, but still, a moment representation similar 



to ( 1.7| was derived in jS80| : suppose the C ^re as above but now two particles coalesce when they have 
spent together an independent exponential time of rate 7. More precisely, suppose w'^ is a solution of 
(1.8) with initial conditions Wq. Then 



i 

where again the product runs over all random walks alive at time t. Sending 7 to infinity for the 
stepping stone model and assuming that Wq = Uq € does not depend on 7, we now observe that 

hm E< [w]{k,) ■ ■ ■ w]{k„,)] = [Utiki) ■ ■ ■ Ut{k^)] , 

7— foo 

since only the coalescence mechanism has changed: random walks now coalesce instantaneously after 
colliding. Boundedness of solutions implies that convergence of the moments suffices to deduce conver- 
gence of the finite dimensional distributions so that the infinite rate limit of the stepping stone 
model is nothing but the standard voter process. For more on this we refer to Section 10.3.1 
of [D9T| . 
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2. Finite Rate Symbiotic Branching Processes 

The aim of this section is to give a compressed overview of definitions and results for symbiotic branch- 
ing processes SBM-y with finite branching rate. After introducing some notation, precise definitions 
and a sketch of existence and uniqueness proofs we turn our focus to the longtime behavior. Let 

(2.1) Ql,:^{W^,W^) 

be the exit law of a pair of (3-correlated Brownian motions started in (w, v) for some u > 0, u > 0, 
stopped at the exit-time 

(2.2) T = M{t:WlW^ ^0}. 

The laws ^ are concentrated on the boundary of the first quadrant which we denote by 

E = {(yi, 0) : yi > O} U {(0, ya) : 2/2 > 0} C 

Whenever the initial condition {u,v) is not crucial we abbreviate the exit-law as Q^. We present in 
the sequel those results on the longtime behavior of symbiotic branching which are related to Q^. 
Those will serve as preparation for the study of infinite rate symbiotic branching processes which we 
wiU denote by SBMoo- 

2.1. Existence, Uniqueness and Tools. Recall Definition |1.7[ where we defined SBM^ as a system 
of coupled stochastic differential equations with drift operator A. With some technical complications, 
Z'' can be replaced by a countable set S and A by an operator 

Aw{i) ^^a{i,j)w{j), 

where (a(i, j)) . ^.^^ is the Q-matrix of a symmetric S"- valued Markov process with uniformly bounded 
jump-rates. The particular case of A occurs for the choice S — 1/^ and a{i,j) = ^ if N ^ j\ = 1- 

2.1.1. State Spaces. Let us define an infinite dimensional state-space for solutions which is commonly 
used in studying interacting particle systems. To do so, suppose /? : — !■ is such that 

^/3(i)<oo and ^ P{i)\a{i, k)\ < M /3{k) 

for all k ^ S. The state-space for the two-type model SBM-y then consists of pairs of sequences that 
grow slowly enough compared to /3: 

Ll = {{u,v) : 5 ^ M+ X M+ s.t. < (X) and < oo}, 

where (/, 5) — J2kes fi^-)9i^)- equipped with the topology induced by the norm ||(u, ■i;)!!^ = 

{\u\ + \v\,(3). Existence of such a sequence f3 is ensured by Lemma IX. 1.6 of |L05| . In the following we 
fix a test-sequence /3 and only work on the corresponding fixed state-space L|. 

2.1.2. Precise Definition and Existence of Solutions. Having defined proper state-spaces, we can give 
the precise definition of solutions to SBM.^. 

Definition 2.1. For {uq,vq) £ L|, we say that {ut,Vt)t>o, more precisely {u,v, , B^) , is a (weak) 
solution 0/ SBM-y on the filtered probability space {^,J^, iJ^t)t>Oi^) 'if 

i) |i3.^(i), B^{i)^ is a set of {Tt)-adapted Brownian motions satisfying for t > 

¥.[Bl{k)Bl{j)\^5^{k-]% 
¥.[B^,{k)B^,{j)\=5,[k~j)t, 
¥.[B}{k)B^,{j)\=Q5o{k-3)t, 
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ii) u.,v. are (J-t)-adapted stochastic processes, almost surely satisfying the integral equations 

ut{k)^uo{k)+ [ Aus{k)ds+ [ y/r^^kyvJk)dBl{k), 
Jo Jo 

Vtik)=vo{k)+ f Avs{k)ds+ f y/jUs{k)vs{k)dB^^{k), 



for k G S, 

iii) {u.,v.) is almost surely continuous with (ut,Vt) £ for all t > 0. 

We now give a quick glance on how to construct solutions for SBM^. For a very detailed proof for 
p = we refer the reader to |DP98| . Prom now on, until the end of the section, we assume 
S — 1/ and ^ = A. The relations to the exit-law of g-correlated Brownian motions remain unchanged 
in this simplified setting so that it serves equally well as a preparation for SBMoo- 

Theorem 2.2. If {uq,vq) E there is a weak solution of SBMj. 

Sketch of Proof. The proof goes along famous arguments due to [SSSOj based on finite dimensional 
SDE theory and limit considerations. Cutting the infinite index set, solutions to the finite system can 
be constructed and then, by moment estimates, their convergence to a weak solution of SBM^ can be 
shown. 

For positive integers n, let Sn [— ti, n]'' be a finite subset of TI^ . To define the approximating 

system, we consider the following system of finite-dimensional stochastic differential equations which 
we denote by SBMi;: 

<(fc) = u^{k) + E ^(""(J') - + /* ^l<{k)vmABY\^), 

<(fc) - v^(k) + f ^(""(J) - ""(^)) ^■'^ + r Vl^^UkH{k)dBf"ik). 

\j-k\=i 



The correlation structure of the Brownian motions remains as in Definition |2.1[ Since this is a system 
of finite-dimensional stochastic differential equations existence of weak solutions 

{unk),vnk),B^'-{k),B'^-{k)}^^^^ 

follows from finite-dimensional diffusion theory for sufficiently "good" coefficients (see for instance 
Theorem 5.3.10 of |EK86j ). To prove non-negativity of solutions, one shows that the semimartingale's 
local time at zero equals to zero (see for instance page 1127 of [DP98j ). 

Solutions (u",u") can be extended to the entire lattice by setting it"(fc) — UQ{k),v"{k) — VQ{k) for 
k Sn- Due to the choice of the initial conditions, the (u",w") are contained in L| for all t >0. 
The main ingredients, to prove convergence of (u", w"), are the following estimates. It suffices to show 
that for fc e 5, T > 0, and e > 

(2.3) supP[sup<(fc) > /sT] -> 0, asK-^oo, 

raeN t<T 

(2.4) sup sup P[K(fc) -<(fc)| > e] ^ 0, as /i ^ 0, 

riGN \t-s\<h,0<t,s<T 



and analogously for u". The desired convergence in ( |2.3| , (2.4) is analogous to (2.9) and (2.10) of 
|SS80j . In order to ensure that all stochastic integrals are martingales we introduce a sequence of 
stopping times: = inf > : -|- > N]. This sequence, almost surely, converges to 

infinity, as N tends to infinity, since solutions do not explode. Using only the definition of SBM" we 
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estimate 



E[ sup 

t<TATV: 



sup 

t<T/\T" 



,13)] -E 

< (uo,/3> +E 



sup 

ktat: 



T - 

^ 2d 

|j-i| = l 



-E 



sup 

t<Tr\T] 



E 



< (?/o,/3) +E 



(2.5) 



TAT" 



E - 

jes„ 

i-ji=i 



-E 



sup /^(*) r \/7<(«)^?(»)dBi'"(*) 



Using the Burkholder-Davis-Gundy inequality and then Fubini's theorem we obtain the following 
upper bound for the above expressions 



E 

»GS'„ 



i-j\=i 



+7E^w 

»GS„ 



So far, this procedure is fairly standard for interacting diffusions of type (1.3 1 where instead of the 
mixed moments, the expectations E[/(u;t(i))] need to be bounded. There, linear growth conditions 
on / lead to a Gronwall inequality which yields the desired bound. In our case, we need to estimate 
moments E[w"(j)] and E[M"(z)D"(i)] uniformly in n. The first moment can be estimated as on page 
1129 of ,DP98, since the correlations do not influence the first moments. The mixed second moment 
is more delicate. Using a point wise representation of solutions, for < 0, the same estimates as in 
|DP98j can be performed. The additional difficulty comes for positively correlated Brownian motions 
{q > 0) that spoil the Gronwall argument in |DP98j due to the appearance of an additional positive 
summand. Nonetheless, the mixed second moment for the approximating system can be estimated 
directly: a moment expression for the finite- dimensional equation, in the same spirit of the moment 
duality that we explain below (see Lemma 2.7), gives the uniform in n upper bound 

(2.6) E[<(*)<(z)] < Ce^^. 

This is similar to the remark on page 41 of |CDG04| where the existence of solutions for g = Q was 
justified by the observation that uv < v? -\- which leads to an upper bound by a two-type Anderson 
model verifying (2.6 1. 



By monotone convergence, getting rid of the stopping times on the lefthand side of (2.5), this implies 

E[sup«,/3)1= lim E[ sup «, /3)] < (uq, /3) + Ct, 

t<T Af-J-oo t<TAT^ 

where Ct is independent of n. Hence, in particular we get by Chebychev's inequality that 

supP[supu"(A:) > K]^Q, 

nGN t<T 



as K tends to infinity To prove (2.4 1 one needs to check that E[|(u" — u", is bounded uniformly 
in n. This can be done similarly as before, using the same bounds on the moments. 
Following the arguments on page 399 of |SS80) . the bounds (2.3), (2.4 1 suffice to ensure convergence 
(in a sufficiently strong sense) of the sequences (it", i;") to a limiting process (u, v) solving the equation 
defining SBM^. □ 
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2.1.3. Tools: Mild Solutions, Total-Mass Processes and Dualities. From the very definition, interact- 
ing diffusion processes are parabolic equations witli random potential functions. In the spirit of the 
deterministic theory one can equally ask for representations that are easier to work with in some 
situations. We will use the weak-solution representation and the variation of constant form. In the 
following we use the semigroup generated by A on Z'^, i.e. the family of linear operators 

(2.7) PJ(fc)= t>0, 

where Pt(i, k) is the transition kernel of a simple random walk on Z'^. 

Convention 2.3. The constant function on taking value u £ is abbreviated by u. 

For a continuum analogue of the following two representations we refer to Corollary 19 of [EF04j and 
for a very detailed proof on the lattice for p = to Theorem 2.2 of |DP98) . 

Proposition 2.4. Suppose that {ut,vt) is a solution of SBM^ with Uq,Vq summable, then Ut and Vt 
are summable and the total-mass processes satisfy 

{ut,l) = {uo,l) + E / VlMj)vs{j)dBl{j), 

(2-8) . 

{vt,l)^{vo,l)+Yl / VlUsU)Vs{j)dBUj), 

where the infinite sums converge in _L^(P). If {uq,vq) G then the point wise representation 



ut{k) = PtUo{k) 

(2.9) 



E / Pt-s{j,kyiUs{])vs{])dBl{]), 



vt{k) ^ PtVo{k) + / pt^s{j,k)^^uS)vs(j)dBl{]), 



holds. The covariation structure of the Brownian motions is as in the definition of SBM-y . 

The weak solution representation can be obtained for {ut,<f)) also for more general test-functions (p. 
Instead of sketching a proof we give an important application leading the way from SBM..y to the 
exit-law defined in (2.1 1. 

A property that is shared by many particle systems is that started at summable initial conditions the 
total-mass process is a martingale. A natural question for the two-type model SBM^ is how the two 
total-mass martingales relate to each other if both types are started at summable initial conditions. In 
order to avoid confusion with (•, •) we denote in the following the cross- variations of square-integrable 
martingales by [•, •]. 

Lemma 2.5. Suppose Uq and Vq are summable, then (wt,!) and (wj,!) are non-negative square- 
integrable martingales with quadratic-variations 

[{u.,l),{u.,l)]^^[{v.,l),{v.,l)]^^j [ {u,,v,)ds 

Jo 

and cross-variation 

Jo 

Sketch of Proof. Positivity of the total-mass processes comes directly from positivity of symbiotic 
branching processes and the martingale property follows from the representation in Proposition |2.4| 
as the martingale property is invariant under L^(P)-convergence. Hence, it suffices to calculate the 
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bracket-processes. The representation stems from the fact that for L'^(P)-convergent martingales also 
the bracket processes converge so that 

[{u.,l),{v.,l)]t= lim [V / ^ju,{k)vMdBl{k), V / ^jUs{k)v,{k)dBUk) 

= Q lim N / 'yUs{k)vs{k) ds 

= £'7 / {us,Vs) ds. 
Jo 

The derivation of the quadratic variations is similar but without the additional correlation parameter 
g. □ 

2.1.4. Dualities. The results on interacting particle systems obtained during the last decades showed 
that the depth of possible results for particular systems depends in many cases on available duality 
relations, i.e. relations of characteristics (here: Laplace transforms or moments) of the process to those 
of other processes. 

In general, one says that a duality between two Markov processes X and Y with state-spaces X and 
y holds if for some duality-function H : X x y ^M. and a potential function f : y ^M. 

(2.10) E''^'[H{Xt,Yo))]^E''"[H{Xo,Yt)efof^''^^'''], t>0. 

In fact, the simplified definition of duality involves potential function f ~ but we will need the 
generalized formulation for one of the dualities for SBM^. Such a semigroup relation can also be 
expressed via generators. Suppose A-^ is the generator for X and the generator for Y, then, under 
some conditions (see Corollary 4.4.13 of [EK86J), (2.10) is equivalent to the generator identity 

(2.11) A''H{;y)^A^H{x,-) + f{-)H{x,-), VxeXyyey. 

The exponential correction is caused by the Feynman-Kac representation for the semigroup with 
additional potential /. The generator relation turns out to be useful in many cases to find a duality 
expression. 

Let us now have a look at the two dual relations for SBM^. For {xi,X2), (2/1,2/2) G L^, define 
(2.12) 

{{xi,x2,yi,y2))e 



E 



- - 0{xi{k) + X2{k)) {yi{k) + 2/2(fc)) + iy/l + g{xi{k) - X2(fc)) (yi(fc) - 2/2(fc)) 



Then, the self-duality relation reads as follows: suppose {ut,Vt) and {ut, Vt) are two solutions of SBM^ 
with initial conditions (uo,fo) € L'p and uq^vq with compact support, then 

(2.13) E"'>^'"'>[exp{{{ut,Vt,uo,io))e)] = E""'"° [exp(((wo, wo, ?it, i^t))^)]- 

For g — the self-duality goes back to |M99j (see also Section 4 of |CDG04| ) and was generalized 
subsequently to £i 7^ in |EF04| . In order to define the infinite rate analogue SBMqo, the self-duality 
will be discussed in more detail in Section 3. On first view the duality looks frightening but it has 
very important applications. Here is an application to weak uniqueness (Proposition 5 of [EF04p of 
SBM^. 

Lemma 2.6. Weak uniqueness holds for solutions to SBM^ if Q ^ (^Ij !)• 

Sketch of Proof. In fact, duality in many cases implies weak uniqueness for corresponding processes: 
for a general result see Proposition 4.4.7 in |EK86| . In this particular case, the proof goes roughly as 
follows: suppose there are two solutions {u\.,vl) and {ul^vD with identical initial condition {uq,vq). 
We aim at showing that both solutions coincide in law. For any fixed pair of compactly supported 
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sequences (j>,ip there is a solution {ut,vt) of SBM^ with initial condition {<f>,tp). Applying the self- 
duality twice shows that for alH > 

E"» [exp {{{ulvlcj,,4,)),)]^E^^^ [exp (((^0,1^0,^^*,^*)),)] 

= E«-''°[exp(((7.?,t;2>,V»,)]. 

A closer look at the complex-valued duality function shows that the equality suffices to deduce the 
claim: in the real direction this is a Laplace transform and in the imaginary direction a Fourier 
transform. Hence, it comes as no surprise that the Laplace transform uniquely determines the law 
of {ul -\- vl, (p + ip) and the Fourier transform uniquely determines the law of (uj — uj, — ip). Now, 
as for finite-dimensional random variables, since (p ip are arbitrary, this uniquely determines 
the laws of the configurations ul + v\ and u\ — v\. Taking sums and differences, the one-dimensional 
marginals of («,*, v^) are determined which finally can be extended to the path- level by Markov process 
arguments. □ 

The second duality is of different type. It does not involve an exponential duality function but a 
polynomial only. The aim is to give an expression for the moments 

(^n+m)] 

in terms of a particle system which we now describe. Suppose that n + m particles in Z"^ are given. 
Each particle moves according to a continuous-time simple random walk independently of all other 
particles. At time 0, n particles of color 1 are located at positions fci, fc„ and m particles of color 2 
are located at positions kn+i, kn+m- For each pair of particles having the same color, one particle 
of the pair changes its color when the time the two particles have spent at same site, exceeds an 
independent exponential time with parameter 7. Let 

— total collision time of all pairs of same colors up to time t, 

hf = total collision time of all pairs of different colors up to time i, 

l\{a) — number of particles of color 1 at site a at time t, 

(a) = number of particles of color 2 at site a at time t 



''t 

and define the duality function 



for (xi,a;2) € i| and Ai,A2 G {a : — >■ N | a(fc) 7^ only finitely many times}. By definition the 
number of particles is finite and constant so that the duality function can be applied to The 
following lemma is taken from Section 3 of |EF04j . 

Lemma 2.7. Suppose {uo,vq) G L| and ki G Z'^, then 
where the dual process behaves as explained above. 



In the general setting A ^ IS., only the dynamics of the single particles needs to be changed: they 
move as continuous-time Markov process with generator A and perform the same changes of colors. 
A simpler situation occurs if the initial conditions are homogeneous. If mq = I'd = 1, then the righthand 
side of the duality only consists of the exponential perturbation. 



Remark 2.8. In the special case of g = 1 and uq = vq — 1, Lemma 2.7 was already stated in |CM94j 



reproved in [GdH07| , and used to analyze the Lyapunov exponents of the parabolic Anderson model. 
Since then the collision times and hf are weighted equally, the colors can be ignored so that only 
exponential moments of collision times need to be analyzed. 
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For p 7^ 1, the difRculty of the dual expression is based on the two interacting stochastic effects: on 
the one hand, one has to deal with collision times of random walks which were analyzed in [GdH07j 
and additionally particles have colors either 1 or 2 which change dynamically. 

2.2. Longtime Behavior and the Exit-Law of g-correlated Brownian Motions. Now, after 
we discussed the basic properties and tools, we can derive the connections of SBM-y and the exit-law 
Q^. This will explain some longtime properties of SBM^. For g = these properties go back to |DP98| 
and [CKOO ■ and they were studied for g G (-1,1) in jBDEllj and for g € (-1,0) in |DM11| . In 
the following we briefly explain the main results on the longtime behavior and sketch the ideas of 



the proofs. Most importantly, we show how the exit law (2.1 1 can be related to limt_>.co (ut, 1) , via 



Lemma 2.5 and then this has several interesting consequences: the weak limit liuit^ca j^iut,vt) can 
be deduced for infinite initial conditions unifying two classical results for the stepping stone model 
and the parabolic Anderson model with Brownian potential. From this, a general technique of jCKOOj 
can be applied to deduce non-convergence in the almost sure sense of the L|-valued processes (ut, Vt). 
As a final consequence of the connection to Q^, we discuss how a critical curve for the behavior of 



higher moments can be deduced which then leads us to the definition of (from Theorem 0.1 1 



2.2.1. Longtime Analysis for the Total Mass - Transience/Recurrence- Dichotomy. A typical question 
for interacting particle systems is the following: does the system get extinct eventually if the initial 
conditions are summable? For several examples this question can be answered easily by taking into 
account nice duality structures. For instance, as explained in the introduction, for the voter process 
this question is equivalent to finite time coalescence of finitely many random walks which again is 
equivalent to recurrence of the random walks. 

The question is more subtle for SBM-^ by lack of knowledge of a useful duality; the self-duality and 
the moment-duality are not very helpful here. Nonetheless, a weaker question can still be answered by 
other arguments. Also note that instead of extinction/non-extinction, the question of coexistence/non- 
coexistence is addressed for two type symbiotic model. Let us first define a notion of "coexistence" in 
the two-type SBM-y model for which we utilize the martingale property of the non-negative total-mass 



processes (recall Lemma 2.5). The martingale convergence theorem implies existence of 



lim {ut, l){vt, 1) =: (itoo, l)(woo, 1) e [0,oo) 
leading to the following definition. 

Definition 2.9. We say that coexistence (of types) is possible, if there are summable initial conditions 
Mo, Wo such that (lioo, l)(i'oo, 1) > with positive probability. Otherwise we say that coexistence is 
impossible. 

For MCBy, Dawson/Pcrkins proved the following recurrence/transience dichotomy. (Recall again that 
we state everything for the case of 5 = Z'' and ^ = A, while some of the results have been proved in 
a more general setting.) 

Theorem 2.10 (Dichotomy for Finite Initial Conditions). Coexistence of types for MCBj is possible 
if and only if a simple random walk on is transient (i.e. d > 3). 

For g ^ the situation is not completely clear, yet. The following result was proved in |BDE11| . 

Proposition 2.11. Let g 0, then coexistence of types for SBM~f[g) is impossible if a simple random 
walk on 1/ is recurrent (i.e. d = 1, 2j. 

For g < the above result has been improved in jPMllj and the following recurrence/transience 
dichotomy has been verified there via moment arguments. 

Proposition 2.12. Let g < 0, then coexistence of types for SBM~f{g) is possible if and only if a simple 
random walk on is transient (i.e. d>3). 
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Propositions |2.11| and |2.12| imply that the question that remains open is whether coexistence of types 
is possible whenever g > and a simple random walk on Z'^ is transient (i.e. d> 3). 

Conjecture 2.13. We conjecture that if d > S and g > 0, there is a critical constant ^{g, d) € (0, oo) 
such that coexistence of types for SBM-y{g) occurs if j < "f{g,d) and coexistence is impossible if 
1 > i{q, d). 

The conjecture is based on a known similar statement for the parabolic Anderson model corresponding 
to £. = 1 (see [GdHOTj). 

In the following we are going to sketch the arguments used in the approach of |DP98| to prove Theorem 

Lemma 2.14. Suppose uq and vq are summable and let T{t) — 7 J^{us,Vs) ds. Then 

{WlWf) ■.= {{uT-^it)A),{vT-Ht)A)) 
is a pair of g- correlated Brownian motions stopped when it hits the boundary E of the first quadrant. 



Sketch of Proof. The claim follows directly from Lemma 2.5 and the Dubins- Schwartz theorem (see 
for instance Theorem 3.4.6 of |KS98| ) applied to both total-mass processes separately. The correlation 
is directly inherited from the driving noises. □ 

With the lemma in hand, let us emphasize the idea behind Theorem |2.10| The pair of total-masses 
is a time-changed planar Brownian motion that stops once it hits E. Hence, in order to prove the 
theorem one has to find a characterization under which the time-change levels off before the 
planar Brownian motion hits E, i.e. one has to show 







Morally, in the transient case islands carrying only type u and islands carrying only type v can move 
away from each other so that {us, Vs) is small even though (wg, 1), {vs, 1) are not. Such scenarios might 
cause the time-change to stop increasing before one of the types vanished. 

Remark 2.15. In fact, the argument of Dawson/Perkins proves more, also for g £ (—1, 1), than .stated 
in Theorem 2.1(\ For d = 1, 2 and g G (—1, 1), the limit ((woo: l)j ("^001 1)) is distributed according to 



the law Q^^y^g 1^ (^y^ 1) of {W^ , W^) (see [BDE11| ) and this will be the crucial ingredient of the proof of 
the next theorem. 

2.2.2. Weak Longtime Analysis for Infinite Initial Conditions in dimensions d = 1,2. An important 
observation of Dawson/Perkins in their study of MCB^ was the use of the self-duality to deduce from 
Theorem |2.10| also the longtime behavior for the system started at infinite initial conditions. This 
task is often much more complicated than in the case of summable initial conditions since simple 
martingale arguments fail. The extension of their approach to 7^ leads to the next result for the 
recurrent regime (see Proposition 2.1 of |BDEllj ). The Brownian motions W^,W'^ and the stopping 



time T are as in (2.1) and (2.2) 



Theorem 2.16 (Infinite Initial Conditions). Suppose that d < 2, g E (—1,1) and uq = u, 
then 

Here, (W^,W^) denotes the pair of functions with constant values {W^,W^). 



Sketch of Proof. Due to the arguments in the uniqueness proof of Corollary |2.6[ it comes as no suprise 
that 

£"^^(ut,vt) converges weakly to C''^''{W^,W^) 
^ lim E"^"[exp(((ut, Vt, (j), ^p))g)] = [cxp{{{W^ , <P, ^)),)], V compactly supported cj), ^. 
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The argument is the same: convergence of the Laplace transform determines weak convergence of the 
sum ut + vt and convergence of the Fourier transform determines convergence of the difference ut — vt- 
Taking sums and differences, convergence of the pair {ut,vt) follows. Employing the self-duality, one 
obtains by dominated convergence 

lim E"'"[exp(((Mt,«t,0,V))e)] 

= lim E^''^[exp(((u,v,7it,^t)),)] 



exp 



( -V^l - giu + v) lim (1, Ut + Vt) + + ~ 1™ (1, Ut - h) 



so that the almost sure conve rgence of the total-mass processes ((1, Ut), (1, vt}) ((1, "oo), (1, Voo}) 
(W^, W^) (see Remark 2.15) gives equality to 



exp 



Finally, it remains to show the identity 



= i?"''^[exp(((VK,i,W^2^0,V')),)] 

for £i-correlated Brownian motions which can either be done via stochastic calculus (for g = see page 
1111 of |DP98| ) or by considering the self-duality for the simplest symbiotic branching model 

\dut = y/'fUtVtdB^ 
l^dvt = y/yurndBt 

(see the proof of Proposition 4.4 of [ BDEll] ). □ 

The need for the assumption d < 2 becomes apparent in the proof: due to the self-duality, the 
convergence claimed in the theorem is equivalent to convergence of laws of the total-mass processes 
to Q^. This occurs precisely in dimensions d— 1,2 (see Remark 2.151. 



The most striking feature of the extension in Theorem 2.16 from g = to g £ (—1,1) is a unifying 
property for SBM-y in the recurrent regime based on the unifying property presented in the Section 
1.2 The theorem is restricted to g € (—1, 1) which is partly caused by the use of the self-duality in the 
proof. However, comparing Theorem |2.16| with well-known results for the boundary cases g = —1, 1, in 
dimensions d = 1,2, it turned out that the classical results can be reformulated in a unified language 
via p-correlated Brownian motions. 



First, suppose wt is a solution of the stepping stone model (see Example 1.4), in dimension d < 2, 
and wq = w £ [0, 1]. It was proved in [S80 that 



(2.14) 



C'^{wt)'^wSi + {l-w)So, 



where Si (resp. 6o) denotes the Dirac distribution concentrated on the constant function 1 (resp. 
0). This can be reformulated in terms of perfectly anti-correlated Brownian motions {W^,W^): For 
g ~ ~1, the pair (14^^, 14^^) takes values only on the straight-line connecting (0, 1) and (1, 0), and stops 
at the boundaries. Hence, the law of {W^, W^) is a mixture of 5(o,i) and ^(i,o) and the probability of 
hitting (1,0) is equal to the probability of a one-dimensional Brownian motion started at w G [0, 1] 
hitting 1 before 0, which is w, and hence it matches (2.14). 



Secondly, let Wt be a solution of the parabolic Anderson model with Brownian potential (see Example 
1.5), in dimension d<2, and constant initial condition wq = w> 0. In |S92] it was shown that 



r'iwt) *^ <5o. 
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As discussed above, if the Anderson model is viewed as a symbiotic branching process with g = 1, 
this imphes 

From the viewpoint of two perfectly positive-correlated Brownian motions in Theorem |2.16| we obtain 
the same result since they simply move on the diagonal dissecting the upper right quadrant until they 
eventually get absorbed at the origin, i.e. {W^,W^) = (0,0) almost surely. 



2.2.3. Almost Sure Longtime Analysis for d 
rem 2.10 states that for d = 1,2 



1,2. Started in homogeneous initial conditions, Theo- 
ut,Vt) converges weakly to a law under which one type completely 
disappears. It is natural to ask whether this convergence also holds pathwise. The negative answer is 
the following result taken from |BDEllj which is based upon the general strategy of jCK00| . 



Theorem 2.17 (Almost Sure Non-Convergence). 
for g e (—1, 1) with uq = u, wq = v. Then 



Suppose (ut, Vt) solves SBM^ in dimensions d = 1,2 



lim inf sup 



ut{k) 
vt{k) 



= 



= 1 



for all {u',v') G E and K C bounded. 



Contrasting the weak convergence to C(W^, W^) of the previous section, the almost sure behavior is 
entirely different: In contrast to choosing, according to the exit-law Q^, one point {u',v') G E as a, 
limit point, the process locally approaches every possible {u',v') E E infinitely often. Hence, looking 
at a fixed box K, the dominant type changes infinitely often and both types approach arbitrarily high 
values. 



Theorem |2.17| is consistent with the known results for the boundary cases even though they appear to 
be very different. Looking inside the proofs of |CKOO| . one realizes that, in fact, they proved more than 
we stated here. Each point of the support of the limit measure (W^, W^) is an accumulation point for 



{ut,vt). Plugging this into the limit measures of the boundary cases discussed below Theorem 2.16 
the theorem extends smoothly to g = —1: The support of u6i + (1 — u)do only contains 1 and leading 



to the classical result that solutions to Example 1.4 alternate locally between and 1 (see Theorem 
2 of [CKOO]). 

Interestingly, the well known almost sure convergence to for the Anderson model from Example 
1.5 (see for instance [GdHOT] ) does not contradict the findings here: the weak limit law is So so that 
the only point of the support is 0. Hence, the techniques of |CKOO| only show that solutions locally 
approach infinitely often which is much weaker than the well known exponentially fast almost sure 
convergence to 0. 



2.2.4. Longtime Analysis for the Moments - The Critical Moment Curve. Here we follow the argu- 
ments of [BDEllj to find bounds for pth moments of SBM-y started at homogeneous initial conditions 
and for pth moments of the total-masses for compactly supported initial conditions. The calculations 
appearing here are the main building blocks for SBMoo so that this section is more elaborate. 

We start with a simple self-duality based lemma connecting the moments of solutions with infinite 
initial conditions with the moments of the total-mass processes for solutions starting at finite initial 
conditions. 

Lemma 2.18. For any keZ"^ andt>0 



(2.15) 



E^'^ [{ut{k) + vtik))P] = E^"'^' [((1, ut) + (1, vt)r] 
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Proof. Employing the self-duality with (f> — ^jj = gives the Laplace transform identity 

where the final equality comes from the uniqueness of the solutions. □ 



With identity (2.15) in hand, to understand the behavior of the moments of the solutions starting 



at homogeneous initial conditions, it suffices to understand the behavior of the moments for the 
total-masses for which we already found a useful structure in Lemma |2.14[ 

A crucial ingredient is an exit-point exit-time equivalence for ^j-correlated Brownian motions (VF^, W'^) 



and their exit time r from the first quadrant Q as defined in (2.1), (2.2). 

Lemma 2.19 (Exit-Point Exit-Time Equivalence). Let p > and u,v > 0, then the following condi- 
tions are equivalent: 

i) 

P < p\q) ■ 



f +arctan(-^)' 
ii) 

£;"'"[r5] < C50, 

iii) 

E''''''[\{W^,W^)\''] < oo. 

The function p is plotted in Figure 1; the critical curve p is strictly decreasing, with 1) — oo, 
p(0) — 2 and p{l) — 1. The particular case of g — corresponds to planar Brownian motion in Q. It 
is a classical result of Feller (see |Fe51| ) that, independently of the initial value, the first hitting-time 
of the boundary only has 1 — e finite moments. This part of the theorem is non-trivial! The second 
part is simpler as the exit point distribution can be calculated explicitly: with 

Q ^ H, z 

the first quadrant is mapped conformally to the upper half plane H C so that, by conformal 
invariance, the planar Brownian path in Q is mapped to a time-changed planar Brownian path in 
H. Luckily, the time-change does not influence the exit-points (only the exit-time) and the exit- 
distribution from H is known to be Cauchy. Plugging this into the conformal mapping, the density of 
the exit-law Q'^ from Q can be calculated (see page 1094 of jDP98) ). The density has no pole at zero 
with tail decreasing polynomially so that the number of finite moments can be deduced. 

For f3-correlated Brownian motions the result follows from a simple change of the space; via 



the £i-correlated Brownian motions are transformed into independent Brownian motions. Simultane- 
ously, the quadrant Q is transformed into a wedge W(0) of angle 



arctan I — , ^ ) 



2 vyr^^ 

using the conformal map 



Z t-^ Z 



The angle of W(0) increases for increasing g explaining, at least morally, the decrease of the number of 
finite moments: if the domain is enlarged, the duration of a planar Brownian path to hit the boundary 
increases, hence, the exit-time might have less finite moments. At the same time, if the Brownian 
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paths run for longer time it will hit larger values so that the hitting-point distribution might have less 
finite moments. Making this rigorous and calculating the exact number of finite moments is done in 
the same manner as for g = 0. That is, the exact distribution of the exit-law can be found: 

P-^'^iW^ = 0, e dr) = ^ , -^dr, 

(2.16) ^ ^li 

P-^^{W^ e dr, = 0) = ^ , 'I ^dr, 



with the constants 

zi = 



n (v — gu)'^\^ /TT / f V — Qu \ ( g 

I -\ — cos — arctan — , +arctan 



(2-17) , , 

zo ^ \ u H — sm — arctan — , + arctan — , 



Q 



From the polynomial decay of the densities given in (2.161 the number of finite moments can be 
deduced. 

Remark 2.20. The explicit density 12.16 ) will play a crucial part in Section 3. For the purposes of 
this section, the density serves as a tool to understand the longtime behavior of moments, whereas for 
SBMryo it will be the main building block of a construction of the process. 

A direct application of the exit-point exit-time equivalence is a proof for the critical moment curve of 
symbiotic branching processes. 

Theorem 2.21 (Critical Moment Curve). Suppose g G (—1, 1) and 7 > 0, then the following hold for 
p > I: 

a) If d = 1,2, then 

P < p{q) *^=^ ^^'^[ut{k)^] is bounded in t > 0. 

b) Ifd> 3, then 

P < p{q) =^ E^'"'' [Mt(fc)P] is bounded in t > 0. 

By symmetry, the same statement holds for u replaced by v. The inverse direction of b) fails and 
depends on 7. 



Sketch of Proof. Taking into account Lemma 2.18 it suffices to prove the equivalence for the total- 
mass process {ut, 1) started at localized initial condition. 

: The proof basically follows from Lemma 2.14 and works for a) and b): The total-masses are time- 
changed £i-correlated Brownian motions and furthermore the quadratic variation (which is nothing but 
the time-change of the Brownian motions) is bounded by r as otherwise one of the total-mass processes 
would become negative. Hence, by the Burkholder-Davis-Gundy inequality 

E^-i'' [((1, Ut) + (1, vt)Y] < CE^'^ . 

The righthand side is independent of t and finite due to the exit-point exit-time equivalence so that 
the claim follows. 

"<^=": Suppose p > p{g). As in the proof of Theorem 2.16 we use the almo st sure convergence of 
{{l,ut), {l,vt)) to {W^,W^) (this only works in case a), see Remark 2.151. Combining this with 
Fatou's Lemma gives 

liminf E^'='^'=r((l,Mt) + (l,Wf))P] > liminfE^"'^" r(l,Mt)P] > E^"^^*- [liminf(l, Ut)^] ^ E^'^\{W^)p]. 
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Figure 1. The critical moment curve p as a function of g ^ (—1,1). 



The righthand side is infinite due to the exit-point exit-time equivalence so that the moment diverges. 
The results for E^'^ [ut{kY] can now be readily deduced by considering the cases ut{k) < vt{k) and 
Vtik) < ut{k). □ 



2.3. Continuuum Model and the Interface Problem. All results discussed above can equally be 
shown for the continuum space analogue model in low dimensions. We will briefly discuss this setting 
as it serves as an important motivation for the study of SBMqo- 

Let us first introduce the model for d = 1. The continuum space symbiotic branching model is defined 
by the pair of stochastic heat equations 



(2.18) 



§iut{x) = 5Aut(x) -I- ^J -fut{x)vt{x) W'^idt.dx), 
§ivt{x) = \/^Vt{x) + ^ jUt {x)vt (x) W"^ (dt, dx) , 



where now A = denotes the typical Laplace operator on M. The driving noises W^,W^ are 
standard Gaussian white noises on x M with correlation parameter g E [—1,1], i.e. the unique 
Gaussian processes with covariance structure 

E[wl{Ai)Wl^{A2)] = (ti A t2)A(Ai n A2), 
E[M/2 (Ai)M/2 (A2)] = (ti A i2)A(Ai n A2), 
E[wl{Ai)W^^{A2)] = g{hAt2)X{AinA2), 

where A denotes Lebesgue measure, Ai,A2 € B{M.) and ti,t2 > 0. Solutions of this model have been 
considered rigorously in the framework of the corresponding martingale problem in Theorem 4 of 
|EF04] ■ which states that, under suitable growth conditions on the initial conditions uo,vq, a solution 
exists for all g G [—1,1]. Uniqueness for g g (—1,1) can be obtained via the self-duality as in the 
proof of Corollary |2.6[ The moment duality also holds with particles moving as Brownian motions 
and collision times replaced by collision local times. 

Stochastic heat equations typically have function- valued solutions only in spatial dimension d = 1. 
The particular symmetric nature of MCB^ changes this property: it was shown in [DEFMPXT] and 
|DFMPX2] that MCB-y do exist in the continuous setting in dimension d = 2 for 7 small enough. 
Existence of solutions in dimensions d > 2 is unknown. 

The results on the longtime behavior will not be repeated here; those are similar to the results discussed 
for the discrete spatial case for d = 1,2. Instead, we include a result of [BDEll refining a Theorem of 
|EF04| . To explain this, the notion of the interface of continuous-space symbiotic branching processes 
is needed. 
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Definition 2.22. The interface at time t of a solution {ut, Vt) of the symbiotic branching model SBM^ 
with p e [—1, 1] is defined as 

Ifci — cl\^x : Ut{x)vt{x) > O}, 
where cl{A) denotes the closure of the set A in M. 

The main question addressed in [EF04j is whether for complementary Heaviside initial conditions 

uo{x) = Ir- {x) and vo(x) = l^+{x). 

the so-called compact interface property holds, that is, whether the interface is compact at each time 
almost surely. This is answered affirmatively in Theorem 6 in |EF04| , together with the assertion that 
the interface propagates with at most linear speed, i.e. for each g £ [—1,1] there exists a constant 
c > and a finite random-time Tq so that almost surely for all T > Tq 

(2.19) U ifct c [-cr,cr]. 

t<T 

For the stochastic heat equation with Wright-Fisher noise corresponding to g = — 1, it was shown 
in |T95j that the correct propagation of the interface is of order VT so that one might ask whether 



(2.19) is sharp for g > —1. Here is a refinement of (2.19), proved in |BDE11) . for which the critical 



moment curve was originally developed. 

Theorem 2.23. Suppose g is chosen sufficiently small such that p{g) > 35 and 7 > 0, then there is 
a constant C > and a finite random-time Tq such that almost surely for all T > Tq. 



U Ifct C 



CVTlog(r), C^TlogiT) 



The strong restriction on g is probably not necessary and is only caused by the technique of the 
proof which is based on the dyadic grid technique utilized for the proof of [T95j . To circumvent the 
boundedness of all moments that holds only for g = —1, 35th moments have to be bounded in time. 
Though the assumption forces g < —0.9958 the result is still interesting. It shows that sub-linear 
speed of propagation is not restricted to situations in which solutions are uniformly bounded as they 
are for g = —1. 

Finally, let us motivate the construction and the study of SBMqo in Section 3. The scaling property 
for symbiotic branching on the continuum (see Lemma 8 of (EF04] ] states that if {ut,vt) is a solution 
started at Heavyside initial conditions, then 

{ut{x),vtix)) := {uNt{VNx),VNt{VNx)) 

is a solution of SBM with Heavyside initial condition. Hence, propagation of the interface of order 

Vt will be intimately related to the behavior of SBM^ with 7 tending to infinity. 
Unfortunately, the constructions in Section 3 can only be seen as a first step towards the correct order 
of interface propagation: the construction for the limiting process SBMqo could so far be carried out 
only for discrete spatial symbiotic branching processes. It is still an open question how to extend the 
characterizations and constructions of SBMqo to the continuum analogue. 
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3. Infinite Rate Symbiotic Branching Processes 



In Section 1.3 we discussed how the standard voter processes can be viewed as an infinite rate stepping 
stone model, or, in other words, SBMoo(f?) f^or g — —1. It is not at all clear if and how that motivation 
extends to g ^ —1 as the coalescing particles duality seems to have no extension to g ^ —1. Taking 
into account the colored particles dual instead, it is by no means clear whether sending 7 to infinity 
leads to a non-trivial process: for 7 = cx) the changes of color occur instantaneously but at the same 
time the exponent is multiplied by cx), so that the moment expression only makes sense if the exponent 
is almost surely non-positive. 

Nonetheless, using the self-duality instead of the moment-duality, it can be shown that sending the 
branching rate to infinity makes sense. To understand the effect in a nutshell, let us take a closer look 
at the non-spatial system of symbiotic branching SDEs 

J dut = ^jutVt dB}, 



(3.1) 



dvt = ^/jutvt dB'^ 



with non-negative initial conditions {u,v). Due to the symmetric structure, we got in Lemma 2.14 
that 

{Wl,W^) := (uT-i(t),«T-Mt)) 

are p-correlated Brownian motions if we use the time-change T{t) = 7 UgVs ds. Caused by the 
product structure of the time-change the boundary E of the first quadrant is absorbing. Hence, the 
Brownian motions W^, W"^ stop at the first hitting-time r of E. Increasing 7 only has the effect that 
(ut, Vt) follows the Brownian paths with different speed so that 7 = 00 corresponds to at once picking 



a point in E according to the exit-measure on E and freeze thereafter (recall (2.1 1). 
To make this argument precise one has to be slightly more careful as the parameter 7 does not 
only occur as multiple in the time-change T but also effects the solution itself. To circumvent this 
obstacle one has to take into account the structure of the equations. Let us label the solutions by 
their fixed branching rate 7. It can be shown that the sequence (u^ ,v^) converges in the so-called 
Me yer-Zheng "pseudo-path" topology (for which we refer to [MZ84] and [J97j) to a limit {U,V). 
Stochastic boundedness in 7 and t of the square-function 7 u2v]ds by t implies that J* UgYs ds = 0. 
Hence, the limiting process {U, V) takes values in E. The only possible limit is the constant process 
{Ut,Vt) — {Uo,Vo),t > 0, where {Uo,Vo) is distributed according to Qg^^ because the prelimiting 
processes {u'^ , v'^) are eventually trapped at £^ at a point distributed according to „. 

Incorporating space, a second effect occurs: both types change their mass on S according to a heatflow. 
This smoothing effect immediately tries to lift a zero coordinate if it was pushed by the exit-measure 
to zero. Interestingly, none of the two effects dominates and a non-trivial limiting process (with 
values in E for each site k S) can be obtained when letting the branching rate tend to infinity. 

Convention 3.1. In contrast to Se ction we do not restrict to the discrete Laplacian A here and 
instead replace A by A as in Section\2.1\ Accordingly, Z'' is replaced by a general countable set S. 



The aim of this section is to explain how the results of |KM11-|-| and |KO10| on the infinite rate 
mutually catalytic branching process MCBgo can be generalized to g 0. After introducing more 
notation for the state-spaces, different approaches to infinite rate symbiotic branching processes are 
presented: a characterization via an abstract martingale problem, two limiting constructions and a 
more hands-on representation via Poissonian integral equations. 

3.0.1. Some Notation. The finite rate symbiotic branching processes were studied on subspaces of 
(M+ X M+) , i.e. at each site of the countable set S the solution processes consist of a pair of non- 
negative values. According to the heuristic reasoning above, at each site k d S infinite rate processes 
take values on the boundary E of the first quadrant so that we can expect to find an £''^-valued 
process. As usual, certain growth restrictions need to be imposed to find a tractable subspace of E^ . 



FINITE AND INFINITE RATE MUTUALLY CATALYTIC BRANCHING 



25 



In accordance with the state-space L| for finite rate symbiotic branching processes we stick to the 
analogue subspace of E^: 

(3.2) lI'^ := Ll n E^ 

equipped with the same norm as L|. Furthermore, we wiU use subspaces of compactly supported and 
summable initial conditions that will be denoted by L^'^ and L^-^ . In contrast to SBM-y, the infinite 
rate processes are not continuous so that solutions have paths in £'([0, oo), L^-^) , the set of functions 
that are right-continuous with limits from the left. 



3.1. Martingale Properties. In order to define infinite rate processes rigorously, in 'KM11-|-| a mar- 
tingale problem characterization was proposed for infinite rate mutually catalytic branching processes. 
This formulation uniquely determines the process but is not very useful for understanding properties 
of the process. Crucial properties of the process, such as non-continuity of sample paths, are not clear 
from this formulation. Nonetheless, it seems to be the most convenient way to introduce the process 
as it directly reveals the connection to the finite rate processes. In what follows we are going to extend 
the results of IKMll-H to Qt^G. 

To define the characterizing martingale problem one crucially uses the self-duality function 
(3.3) F{xi,X2,yi,y2) = cxp {{{xi,X2,yi,y2))e) 



defined in (2.12). We include the next two simple (stochastic) calculus lemmas in order to clarify the 



appearance of g in the definition of (( 

-2 

-F{xi,x2,yi,y2) = F{xi,x2,yi,y2){{l,0,yi,y2)) s 



Lemma 3.2. Suppose {xi,X2) G L'i and (2/1,^2) compactly supported, then for all k G S 



d 

dx\{k) 
d 

dx2{k) 



F{xi,X2,yi,y2) = F(a;i, 12, y2)((0, 1, yi, 1/2))^ 



and 



1 

2dxi{ky 



1 92 



92 



2dx2{kY ^ dxi{k)dx2{k) 



F{xi,x2,yi,y2) = 4(1 - g'^)F{xi,x2,yi,y2)yi{k)y2{k), 



where d/dxi{k) (resp. d/dx2{k)) denotes the partial derivative with respect to the kth coordinate of 
the first (resp. second) entry. 

Proof. First note that all appearing infinite sums are actually finite as yi and j/2 are compactly 
supported. We leave the simple derivations of the first derivatives to the reader as it does not clarify 
the influence of g. 

Abbreviating c(fc) — yi{k) + 2/2 (^) and d{k) — yi{k) — y2{k), by the chain rule we obtain 

^1 92 ^ Q2 92 " 



2dxi{kY 2dx2{kY dxi{k)dx2{k) 



F{xi,x2,yi,y2) 



F{xi,x2,yi,y2) 



^ ( - v/r^c(fc) + i^/lT~gdik)^ ' + K " yi^c(fc) - iyiT~gdik)^ ' 



+ g[ - - gc{k) + Vl + gd{k)j ( - ^/l - gc{k) - i^\ + gd<yk)^ 
which is equal to 

F{x^,X2. yi, ya) [(1 - g){^ + g)c^{k) - (1 -f g){\ - g)d^{k)\ 
= -F(xi, 0:2, 2/1, 2/2)4(1 - ^)yi{k)y2(k). 



□ 
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The intrinsic need for the particular choice of ((•, •, •, ■))g can now be revealed: the additional square- 
roots involving g are chosen in such a way that the cross-variations caused by the correlated driving 
noises cancel. 

Proposition 3.3. Suppose (uo,fo) G (2/1,2/2) G L^'^ and {ut,Vt)t>o a symbiotic branching 
process with finite branching rate 7 > and correlation parameter p G [—1, 1], then 

(3.4) ft 

F{ut,vt,yi,y2) ~ F{uo,vo,yi,y2) - / {{Aus,Avs,yi,y2))gF{us,Vs,yi,y2)ds, t>0, 

Jo 

is a martingale null at zero. 

Proof. Noting again that all infinite sums are, in fact, finite as the test-sequences yi and j/2 have 
compact support, we may apply Ito's formula to the finite set of stochastic differential equations to 
get 

F{ut,vt,yi,y2) 
^ F{uo,vo,yi,y2) + / ^ — 77TF{us,Vs,yi,y2) dus{k) 



+ ^E/ d£(ky2Fius,Vs,yuy2)d{v.{k)), + Y,J^ ^("^. ^i. y2)d(^..(fc), ..(fc))., 

where we used that by definition the Brownian motions at different sites are independent. The corre- 
lation structure for the Brownian motions at the same sites and the previous lemma yield equality of 
the above expression to 

F(uo, Wo, 2/1,2/2) -I- local mart. + / F{us, Vs, yi, y2){{Aus, 0, yi, y2)) g ds + local mart. 

Jo 

+ / F{us,Vs,yi,y2){{0,Avs,yi,y2))gds + / F{us,Vs,yi,y2)^^il - g'^)yiik)y2ik)^Us{k)vs{k) ds. 
Sorting the terms leads to 

F{ut,vt,yi,y2) = ^"("0, ^o, 2/i, 2/2) + local mart. + / F{us,Vs,yi,y2){{Aus,Avs,yi,y2))eds 

Jo 

+ / F{us,Vs,yi,y2)^^il - g'^)yiik)y2{k)^Us{k)vsik) ds. 
■^0 fees 

By assumption y e L-f'^ so that yi{k)y2{k) = for all k & S. Hence, the last summand vanishes 
and it only remains to show that the local martingale is a martingale. But this follows directly from 
the fact that F is bounded and the moments E[ut{k)vt{k)] are locally bounded. The latter follows for 
instance from the moment duality of Lemma |2.7[ □ 



It would be desirable to uniquely define solutions of finite rate symbiotic branching processes via this 
martingale property which unfortunately is impossible: the corresponding martingale problem does 
not involve 7 and it is satisfied by SBM-y(£i) for arbitrary 7. As symbiotic branching processes for 
different branching rates do not coincide in law, the martingale problem has infinitely many solutions. 
However, the class of processes on the restricted state-space i^'^ is less rich so that the small class of 
test-functions suffices here for the martingale problem to be well-posed. In particular, the restriction 
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rules out all solutions of SBM^. Here is the generalization from g = to g E (—1, 1) of Proposition 
4.1 of IKMII+I . 

Theorem 3.4. Let g e (—1, 1), then there is a unique solution to the following martingale problem: 
For all initial conditions {xi,X2) € L^' , there exists a process {U,V) with paths in Z?([0, oo), L^' ) 
such that for all test-sequences (yi,2/2) G L^'^ the process 

(3.5) Mt°°{x,y) :^ F{Ut,Vt,yi,y2) - F{xux2,yi,y2) - f {{AUs,AVs,yi,y2))eF{U,,Vs,yi,y2)ds 

is a martingale null at zero. The induced law on ^([O, oo), L^'^) constitutes a strong Markov family 
and the corresponding strong Markov process will be called infinite rate symbiotic branching SBMao{g)- 



We postpone a sketch of a proof to Section 3.4 where solutions are constructed by means of the 



Poissonian equations already mentioned in Theorem 0.1 

Since we discussed extensively the longtime behavior of finite rate symbiotic branching processes we 
say a few words about the longtime behavior of infinite rate symbiotic branching processes. The case 
oi g — Q has been studied in |KMfObj and some sufficient conditions for coexistence and impossibility 
of coexistence have been derived there. For p < a full recurrence/ transience dichotomy has been 



established in |DM11) in the spirit of the results presented in Section 2.2.1 



Proposition 3.5. Let g <Q, then coexistence of types for SBMaoig) is possible if and only if a Markov 
process on S with Q -matrix A is transient. 

Note that this proposition extends Proposition |2 . 1 2 1 to 7 = oo on a general countable site space S and 
an arbitrary symmetric Markov process with Q-matrix A. For the proof we refer the reader to jPMll] . 

3.2. Main Limit Theorem. So far we have discussed the finite rate symbiotic branching processes 
and introduced the well-posed martingale problem from which one can define the family of processes 
SBMoo(£'), g € (—1, 1). To get the link between the two, we sketch in this section how to show that 



SBM.y converges in some weak sense to the solution SBMqo of the martingale problem (3.5 1 as 7 goes 
to infinity. This, in fact, justifies to call the processes of Theorem [T4] infinite rate symbiotic branching 
processes. 

Unfortunately, the convergence of SBM.^ to SBMoo will not hold in the convenient Skorohod topology in 
which continuous processes converge to continuous processes. As a solution of the system of Brownian 
equations (0.1 1, SBM^ is continuous, whereas SBMqo is non-continuous as solution to the system of 



Poissonian equations. 

Even though the convergence can not hold in the Skorohod topology, it holds in some weaker sense. The 
suitable "pseudo-path" topology on the Skorohod space of RCLL functions was introduced in |MZ84] . 
The topology is much weaker than the Skorohod topology and is, in fact, equivalent to convergence 
in measure (see Lemma 1 of |MZ84| and also results in (J97) ). Sufficient (but not necessary) tightness 
conditions for this "pseudo-path" topology were given in |MZ84j . In particular, these conditions are 
convenient to check the tightness of semimartingales. 

Here is the extension of Theorem 1.5 of |KMlf+1 to 7^ 0. 



Theorem 3.6. Fix any g G ( — 1, 1). Suppose that for any 7 > 0, {u] ,Vf)t>o solves SBMy{g) and the 



initial conditions (uq,wJ) = (C/o,Vo) G do not depend on 7. Then, for any sequence 7„ tending 

to infinity, we have the convergence in law 

{u'^^v'^^') =^ {U,V), n^oo, 

in D{[0, 00), L'p) equipped with the Meyer-Zheng "pseudo-path" topology. Here, (U,V) is the unique 
solution of the martingale problem of Theorem\3.4\ 
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Sketch of Proof. The proof consists of three steps: 

Step 1: Tightness in the Meyer-Zheng "pseudo-path" topology follows from the tightness criteria of 
|MZ84| . To carry this out, one has to show tightness for the drift and the martingale terms in the 
definition of SBM-yi By standard estimates the drift terms are, in fact, tight in the stronger Skorohod 
topology: this follows from 

(3.6) supE[sup(u7,/3)H < oo, pe{l,p{g)). 

Apart from the facts that /3 ^ 1 and Uq is not assumed to be summable this is close to the moment 
bounds for the total-mass processes that we obtain from Lemma |2.18| and Theorem |2.21| With the 
same trick as in Lemma 6.1 of |KM11-|-| , the lefthand side of (|3.6[ ) can be bounded uniformly in 7 by 
a multiple of E[t^'/^], where t is the exit-time of Theorem |2.21[ Replacing p € (1, 2) by p e {l,p{g)), 
the arguments in the proof of Lemma 6.2 of iKMll-|-| carry over line by line. The crucial observation 
is that for all g £ (—1, 1) the critical curve p{g) is strictly larger than 1 which is all that is needed. 
To prove tightness of the martingale part, the tightness criteria for martingales can be applied (compare 
Theorem 4 combined with Remark 2 of ^MZ84) ). 

Step 2: To show that all limit points indeed solve the martingale problem, we only have to use 



Proposition 3.3 and some moment estimates. For all fixed 7 > the same martingale problem is 
fulfilled so that it comes as no surprise that the martingale problem is fulfilled in the limit if one can 
show that the martingales converge to a martingale. But this follows from the same estimates that 
are used for the tightness proof involving crucially the critical curve p{g). 



Step 3: In the previous section we stressed out that the martingale problem (3.5) is only well-posed if 



the involved process takes values in the restricted space L^p^ . To show that for any limit point {U, V), 
we indeed have Ut{k)Vt{k) = for alH > 0, A: e S, one can show that almost surely 

(3.7) / {Us{k)Vs{k) Al)ds ^0, Vt > 0, fc e 5, 

Jo 

since then, by right-continuity, {Ut{k),Vt{k)) G E for all i > 0, fc G S. By tightness of step 1 one 
can easily derive stochastic boundedness of 7 /p ujv] ds uniformly in 7 for any A: g S" as in the proof 
of Lemma 6.3 in t KMll-|-| from which ( |3.7| follows by taking into account that convergence in the 
Meyer-Zheng "pseudo-path" topology is equivalent to convergence in (Lebesgue) measure. □ 

Now that above we have made precise sense of SBMqo in terms of a weak limit of SBM^ that solves 
a well-posed martingale problem, the next two sections are devoted to constructions that shed more 
light on the properties of the processes. 

3.3. Trotter Type Construction. A very different perspective for MCBoo was presented in |KO10| . 
Their main idea was to combine "by hands" the precise infinite rate limit for the mutually catalytic 



SDE (3.11 with the heatflow corresponding to the generator A, to construct a more instructive ap- 
proximation. The approximation converges in the stronger Skorohod topology, instead of only in the 
weaker Meyer-Zheng "pseudo-path" topology, which might be helpful to deduce properties for the 
limiting process. We now briefly discuss here how their approach extends to SBMoo(£') for g E (—1, 1). 
Separating the deterministic and stochastic terms in the very definition of SBM.y, one has to consider 
the pair of evolution equations 

d d 

(3.8) a^"* ^ Wt^^^ 

and the set of independent two-dimensional symbiotic branching processes 

(3.9) dut{i) = \/'^ut{i)vt{i) dBl {{), dvt{i) = jut{i)vt{i) dBf (i) . 



The evolution equations (3.8) can be solved explicitly in terms of the semigroup Pt corresponding to 



A (recall (2.7 1 for A = A) and the solutions do not depend on the branching rate 7. The processes 



in (3.9) obey a more interesting behavior as we have discussed in the introduction of this section: the 
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pairs of independent stochastic integrals provide a set of independent diffusions indexed by S which, 
for 7 = oo, correspond to a set of independent choices of the exit-law Q^. 

The Trotter type approach to SBMqo is built upon these two explicit representations: the heatflow 
based on A and the independent choices of the exit-law are alternated with increasing frequency. 
For e > 0, the approximating processes {U'^,V'^) are defined as follows: 



i) Within each interval [ne,{n+ l)e), {U^,V^) is the explicit (deterministic) solution of (3.81 
with initial condition (U^^, V^^). 

ii) At times ne, {U^^_, V^^_) is replaced at each site k € S independently by a point in E chosen 
from the exit-law of p-correlated Brownian motions started at (U^^_{k),V^^_{k)) . 

It follows almost directly from the definition of the approximation that for any e > fixed, {If^, V^) 



is a solution to the martingale problem (|3.5|). As discussed below Proposition 3.3 this does not cause 

-2, J 



any contradiction since even with initial condition {Uq, Vq) G L'I'^ , the processes {If^jV^) take values 



in L| but not in the restricted state-space L^' . The adaption of the main result of |KO10| is then as 
follows. 

Theorem 3.7. Suppose g G ( — 1,1) and {Uo,Vo) G L^p^ ■ Then, as e tends to zero, the family 
{{U'^,V'^)}e>o converges weakly, in the Skorohod topology on £>([0, cx)), L|), to SBMoo{g). 

Sketch of Proof. The proof consists of three steps: 

Step 1: Tightness in the Skorohod topology is proved via Aldous' criterion and moment estimates 
that are based on estimates for the exit-measures Q^. The extension from g — to g € (—1,1) is 
crucially built upon the fact that the estimates of [KOIOJ are based on boundedness of some moments 



greater than 1 and this equally holds for any g € (—1, 1) (see Lemma 2.191. 

Step 2: The identification of the limit points is not difficult since the approximating sequence already 
solves the martingale problem for any e > 0. It is only needed to show that the sequence of martingales 



remains a martingale for which again moment estimates based on Lemma 2.19 are needed. 



Step 3: The limiting process takes values in the smaller space L^^'^ C due to the construction and 
continuity of the heatflow. □ 

3.4. Poissonian Construction. Up to now the infinite rate symbiotic branching processes have only 
been characterized as weak limits of approximating sequences and via an abstract martingale problem. 
The most explicit construction of SBMqo is presented here as the unique weak solution to a system 
of Poissonian integral equations from which we deduce the connection to the voter process. 



3.4.1. Jump Measure. To describe the jumps of SBMqo, the following definition is needed: 

Definition 3.8. Suppose Q^ u,v) ^'^ exit-measure of g- correlated Brownian motions started at {u,v) 
from the first quadrant (see (2.1)). Define 



where the limit is in the vague topology on measures ( i. e. integrated against continuous functions with 
compact support). 

Recalling that by definition is a probability measure for arbitrary initial conditions, after rescaling 
the measure has to be an infinite measure on E equipped with the restricted Borel cr-algebra. 

Next, a density for will be derived. The core of the work has already been done in |BDE11| where 
the density for was calculated (see (2.16)); with this density in hand we get the following result. 
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Lemma 3.9. For g e (—1,1) and a > the measure is absolutely continuous with respect to the 
two-dimensional Lebesgue measure restricted to E with the density 



^(a,0)(%l'y2)) 

where p(g) — and 9(g) — ^ + arctan 



dyi ■y2 = 0, 
dy2 : 2/1 = 0, 



(x/i^)- 



Proof. By definition of ^ all we need to do is to plug-in = (i, e) into the the explicit density 

given in (2.16), divide by e and go to the limit. With the notation used in (2.16), (2.171 we obtain 



( arctan 



ag 



g^a 



+ arctan ( 



Q 



e— > 

a- 







p(q) 



Taking into account sin(2:) ~ a; at and I'Hopital's rule with arctan'(a;) = 1/(1 + x^), leads to 

^2 



{e-agy\2wt^ . ( n 

— I sm I 

1 - g^ 



sm I —— i arctan 
\9{g)\ 



ag 



-arctan ( 



9 



Vl-< 



2«(e) _i . 

e sm 



p(e) TT 



0(e) 



Plugging this calculation into (2.16), (2.17) the claim follows. 



□ 



In fact, it will always be sufhcient to consider the case a = 1 by simple scaling as we will see below in 
Lemma 13.111 

Definition 3.10. The special case for a — 1 will serve as basic jump measure. We abbreviate 

„p(e)-i 



v%d{yi,y2)) = i^^^Q^{d{yi,y2)) = < 



P(g)^v/r^ /^l , dyi :y2=0, 

7r(yi "'-Ij 

p{q?VT^ /ilt\'^ dy2 :yi=0. 



in the sequel. 



A quick glance at the explicit density of i^^ shows that the densities are far from being symmetric for 
the X- and y-axis: a pole is found only at (1, 0). Moreover, the tail behavior shows that the measure 
restricted to the y-axis is finite and is infinite for the x-axis. This comes as no surprise from the def- 
inition: starting the correlated Brownian motions in (l,e) and sending e to zero forces the Brownian 
motions to exit the first quadrant closer and closer to the point (1,0). The additional factor 1/e then 
leads to the pole at (1, 0). 



The reduction from v 



ui^ 1^1^-^ is motivated by the following scaling property. 



to l'^ 

Lemma 3.11. Suppose f maps E to R continuously, then 
I f{yi.y2)vtaM{yi,y2)) = -[ 

J E 0. JE 



f{ayi,ay2) {d{yi,y2)). 
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Proof. Splitting E in the two positive parts of the axes, the claim follows from a change of variables 
in the third line of the following computation: 



/(yi,0);/(%)(%i,0))+ / /(0,y2)<,o)('^(0'y2)) 

E J E 

, ,,p(e)-i 

- /(yl,0)p(^.)2v/^^a^(^)-l— |i —^dyr 

T^{yi - aP(e)) 

Jo Tr{y^^^' + aP(e}^ 



,,p(e)-i -I /-oo p(e)-i 

"2 



1 

a Je 



f{ayi,ay2)i^'^{d{yi,y2))- 



□ 



3.4.2. Poissonian Integral Equations. The aim of this section is to discuss properly the objects ap- 
pearing in Theorem |0.1| and to give elements of the proof. For convenience of the reader not familiar 
with jump diffusions we added a (very brief) summary to the appendix. 
Let A/" be a Poisson point process on S x E x (0, oo) x (0, oo) with intensity measure 

(3.10) Af\{k},diy,,y2),dr,ds) = i^<^{diyuy2))drds, Vfc e 5. 

The Poisson random measure Af can be interpreted as a collection of independent Poisson point 
measures {Af{{k},-,-,-),k G S} on E x (0,oo) x (0, oo) running independently at each site k G S. 
Then at each site k Cz S, the basic jump measure will be used to determine the target point of 
a jump from _E to at that site. To incorporate a state-dependent jump rate the r-component will 
be used. That is, as the jump intensity will depend on the current state of the system before time s, 
which we denote by {Us-,Vs-), we define the intensities 



(3.11) Is{k) 



■Us-{k)>0, 

v;-(fc) >o, 



V,-{k) 





Us-{k)^Vs-{k) 



Let us take a closer look at (3.111, and assume for a moment that the current state at site k is 
(Us-{k),0). Then the intensity of jumps at k at time s is high if Us-{k) is small compared to the total 
size of the population of "type y at neighboring sites. 

Next, we need to specify the integrand that describes the jumps of SBMqo at an atom of M at 
{k, {yi,y2),s,r): 
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SO that at an atom (fc, {yi, 1/2), s, r) of J\f the system in state {Us-, Vs-) changes at site k via one of 
the foUowing transitions: 



(3.12) 



(k) 



i ' 

\y2Us- 

( ° 

\yiVs- 



y 



('.■) 



y = 



The second and fourth cases will be referred to as change of type as the jump changes the current 
state from one axis to the other. Next, after compensating the jump term and adding an additional 
drift term, we are ready to define the system of Poissonian equations: 

Definition 3.12. The system of Poissonian integral equations, indexed by the possibly infinite set S, 

ft 

AUs{k) ds 



Ut{k) = C/o(fc) 

r-i i-I^ik) 







(3.13) 



/o Jo 

Vt{k) = Vo{k) 



hik) 



y2Vs-{k) + (2/1 - !)[/,_ (fc))(AA - AA')({fc}, d{y^,y2).dT, ds), 
AV,{k) ds 

y2Us-{k) + (yi - 1)K- (fc)) (AA - ^f')m. diyi,y2),dr, ds). 



in short 

Ut{k) 
Vt{k) 




J^AUsik)ds\ 
J^AVs{k)ds) 

Vs-{k) 
Vs-ik) 



2/2 



+ (yi - 1) 



Us-{k) 

Vs-{k) 



{M -M'){{k},d{yi,y2),dr,ds) 



will be called infinite rate symbiotic branching SPDE. 



For a further discussion and connections of the above Poissonian SPDE to the standard voter process 
we refer to the next section. 



Theorem 3.13. Let g e (—1,1) and suppose {Uq,Vq) £ L^p^ . Then Equation (3.13) admits a weak 
solution with paths almost surely in Z)([0, cx)), L^'"^). 

Sketch of Proof. The proof is along the lines of Sections 2 and 3 of |KM11+| for gi = 0. The basic idea 
is to construct the approximating sequence of equations with the following modifications: 

(1) truncate the infinite index set (compare with the proof of Theorem 2.2)), 

(2) modify the jump measure by truncating its jumps near the pole (1,0). This makes the 
modified jump measure finite. 

(3) modify the jump intensity / by truncation its big values. 

With these truncations there are only finitely many jumps up to any t > so that solutions can be 
built merely "by hands" via interlacing. To be more precise, we consider equations on subsets S"" C S 
with m elements and we redefine 



v^'^{d{yi,y2)) = v^{d{yi,y2))l{y^-i>,.i^y^>e'}- 
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This asymmetric truncation around (1,0) is slightly strange but if e, e' are chosen such that J^{yi — 
l)i''^'^{d{yi,y2)) = then solutions stay on the boundary E of the first quadrant, i.e. the drift does 
not push solutions into the interior. The modified state-dependent jump rate becomes 



Itik) = 



C/j"l''(fc)Ve 
-4C/"''(fc) 
V™''(fc)Ve 

\AUrj\k) + \AVll-\k) 



(k) > 0, 
(fc) > 0, 



ur-\k) = v,"i^\k)^o. 



Replacing ly^ by i^'^'^, I by P and adding additional jumps away from (0, 0), solutions ([/'"''^, V"^''^) can 
be constructed by hands via a Poisson point measure Af since jumps do not accumulate. By definition 
it seems clear (ignoring the additional jumps away from zero which will vanish in the limit) that a 
possible limit for m — > oo and e — > fulfills Equation (3.13). This can be made rigorous via the method 



of characteristics for semimartingales and classical convergence theorems. To ensure that the sequence 
converges, tightness in the Skorohod space is justified by Aldous' criterion. Up to now the arguments 
copied directly those of IKMll+l , with only difference in replacing i^^ by i^^. To apply Aldous' criterion 
one then has to find pth moment estimates which again can be obtained as in IKMll+l by replacing 
their arguments relying on p e (1, 2) = (l,p(0)) by the same arguments based on p G {l,p{g)). □ 



In order to clarify the connection of Equation (3.131 to infinite rate symbiotic branching processes 



we give more elaborate arguments for the next result which for £i = was proved in Lemma 3.12 of 
IKMll+l . 



Proposition 3.14. Let g G (—1, 1) and suppose {U, V) is a weak solution to Equation \3.13 ) taking 
values in L^^'^ , then {U, V) is a solution of the martingale problem (3.5). 



Sketch of Proof. To show that the martingale problem (3.5) is satisfied by weak solutions to (3.13), 
one can proceed similarly to the proof of Proposition 3.3 by applying Ito's formula to i^(J7t, Vj, zi, 22) 



for compactly supported zi,Z2. Let us first define the integrands of the Poissonian integrals as 



' .h{yi,y2.Us-{k).Vs-{k)y 
J2{yi,y2,Us-ik),V,^ik)) 



2/2 



Vs-{k) 

Us-{k) 



{yi 



Us-{k) 

Vs-{k) 



and abbreviate for (a; 1, 0:2), (^i, 2^2) G 



{{xi,X2,Zi,Z2))g^k 

= -v/1 - Q{xi{k) + X2{k)) (zi(fc) + Z2(fc)) + i^/l + Q{xi{k) - X2{k)) (zi(fc) - Z2{k)). 
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First, by Ito's formula for non-continuous semimartingales and the notation for partial derivatives 
already used in Lemma 3.2 we obtain for {zi, Z2) G L^'^ 



F{Ut,Vt,zi,Z2) 
= F{Uo, Vo, z,,Z2) + E /* ^^/'^''-''-''^''^^^AUM ds 



E 

kes 

E 



d 



fees 



dx2ik) 



^0 



MiU..Vs) + J{Vl.y2,U,.y,),Zi,Z2))i,,k _ p{{U,,V,,ZuZ2)),,k 



Jiiyi,y2,Usik),Vsm 



J2iyi,y2,Usik),Vsik)) 



dU{k) 



dV{k) 



A{Us,Vs,Zl,Z2))i,,k 



MUs.Vs,Zi,Z2))s,k 



J\f'{{k},d{yi,y2),dr, ds) 



local martingale 



which, carrying out the partial derivatives via Lemma 3.2 and plugging-in the definition of A/"', yields 
F{Ut,Vt,zi,Z2) 

= F{Uo,Vq,zi,Z2) + / F{Us,Vs, zi, Z2){{AUs, AVs, zi, Z2)) g ds + local martingale 
Jo 



- {{J{yi,y2,Us,Vs),zi,z2))e,k 



■,{{J{yi,y2,Us,Vs),zi,z2))g,k 



v^{d{yi,y2))ds. 



The righthand side is already close to the martingale problem (3.5 1 if we can show that the sum of 
the integrals with respect to v^{d{yi^y2)) ds equals to zero and the local martingale is a martingale. 
Note that this is very similar to the proof of Proposition |3.3| To prove the first assertion, note that, 
by definition of , 



e^Wy^^y-'^-''^^'^^^^-^^^^'' ^1- {{J{y,,y2,Us,V,),z,,Z2))g,k v'{d{y^.y2)) 

^ " [e^^O^.^^.^-^^.v-^)^--^^))..^ - 1 _ {{J{w:^,W?,Us,Vs),z,,Z2))g,k 
so that we are done if we can show that, for any (xi, X2), {zi, Z2) S and e > 0, 



lim 



(3.14) 



= 0. 



But this identity holds, if r is replaced by < > 0, by Ito's lemma as in the proof of Proposition |3.3[ 
The necessary arguments that justify the changes of limits and integration are as in Lemma 3.11 
of IKMll+l . Those incorporate the exit-time exit-point equivalence of Lemma 2.19 for p 7^ 0. The 
martingale property for the local martingale then follows as in the proof of Lemma 3.12 of IKMll-H 
from first moment estimates that are not affected by € 1)- D 



In the proof Proposition 3.14 we did not utilize the particular form of the intensity /t(fc) so that 
arbitrary changes in the intensity seem to lead to other solutions of the martingale problem, and by 
this seemingly imply a contradiction to uniqueness of the martingale problem (3.5 1. However, this 
chain of reasoning is not true because of the particular choice ( 3.11 ) for /t(fc) forces solutions to have 
paths in L^' and uniqueness for the martingale problem (3.5) only holds for solutions with paths 
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restricted to i^'^- 

Let us make this more precise: suppose that Uxik) — for some random time T > and some fc e S*. 
From the density of the basic jump measure i'^ it is clear that for some positive time no jump changing 
the types occurs (by finiteness of restricted to the y2, jumps that change types come with finite 
rate). Hence, for some positive random time S, no jump occurs so that 

C/T+r(fc) = 0, re [0,(5]. 

In particular, this shows that Af' must be such that 

rT+r pT+r ph{k) p 

AUs{k)ds^ / / Ji(2/i,y2,f/.-(A:),-K,-(fc))AA'({fc},d(yi,2/2),dr-,ds) =0, 



for all r £ [0, S]. We now briefly show that the choice (3.11 ) indeed does the job: 

rT+r 

' Ji{yi,y2,Us-{k),Vs-{k))I,{k),y<^{d(yi,y2))ds 



T J E 

T+r 



T Jo 

T+r 

AU.^ik) ds 



All (k) , 



dy2 ds 



because 



y2P 



i/2 



„p(q) 



dy2 = lim -E^'^\W^] = lim - lim E'^■■'\W^.^ = lim -e = 1. 



Note that here the superscript in W refers to the second coordinate of the pair of Brownian motions 
and not to the second moment. The first equality follows from the definition of the second follows 
from the martingale convergence theorem for which the uniform integrability is ensured by the upper 
bound 

where the positive constant ^ is chosen sufficiently small such that p{q) — > 1 (existence of fx is 
ensured by the exit-time exit-point equivalence of Lemma 2.19). 



With the Poissonian construction of SBMqo in hand we now sketch a proof of Theorem 3.4 



Sketch of Proof for Theorem \3.4\ Existence of solutions to the martingale problem follows from The- 
orem |XT3] and Proposition |3.14 



The uniqueness proof is inspired by the proof of Lemma 2.6 for 7 < oo based on self-duality. Here, 
we sketch the chain of arguments of Section 4 in |KM11-|-| which can be copied line by line while 
replacing the duality function in |KM11-|-| by the p-dependent duality function F defined in (3.3). 
Step 1: For compactly supported initial conditions (C/g, Vq) solutions (f/, V) to the martingale prob- 
lem are constructed via the Poissonian equations ( 3.13[ ). From the first moment estimates one obtains 
that solutions decay sufficiently fast at infinity. 

Step 2: First moment bounds for arbitrary solutions of the martingale problem are derived by differ- 
entiating the Laplace transform part (see Lemma 4.2 of ,KM11+. for g = 0). 
Step 3: The crucial part is to derive the self-duality relation 

nF{UuVu ilo, Vo)] = E[i^(t/o, Vo, UuV)] 



between the two independent solutions {U,V) and {U,V) starting at {Uo,Vo) € L'i'^ and {Uq,Vo) E 



L^'^ . Now, as in the proof of Corollary 2.6 self-duality determines the one-dimensional laws {Ut, Vt) 
along the lines of the proof of Proposition 4.7 in [KMll-l-J for g — 0. Standard theory (see Theorem 
4.4.2 of |EK86j ) allows us to extend the uniqueness of 1-dimensional distributions to uniqueness of finite 
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dimensional distributions. Finally, the strong Markov property for (C/, V) follows from measurability 
in the initial condition which is inherited from the finite jump rate approximation processes. □ 



Combining Theorems 3.4 3.13 and Proposition 3.14 we immediately get the following theorem. 



Theorem 3.15. Let g G (—1, 1) and {Uq, Vq) G L^p^ ■ Then there exists unique weak solution to (3.13) 



which is the unique solution to the martingale problem from Theorem \3.Ji\ 



3.5. Infinite Rate Symbiotic Branching Processes and Voter Processes II. The infinite 
rate symbiotic branching processes SBMoo were characterized in previous subsections via various 
approaches. In this final section we describe SBMoo from the viewpoint of the standard voter process 
which is closely related to symbiotic branching with g = — 1 as we have already seen in the Section 

For the rest of this section we stick to ^ = A on S* = Z'' for convenience. 



We start with restating Theorem 3.6 for the g — —1 case. However, note that we additionally have 
to assume Uq + = 1 since we cannot use the self-duality anymore as for g — —1 it does not carry 
enough information to characterize the full law of the limiting process {U,V). Under this additional 
assumption we can rely on the folklore results mentioned at the very end of Section |1.3| whereas for 
general initial conditions a different approach should be developed. 

Theorem 3.16. Suppose g — —1 and for any 7 > 0, [u] ,v])t>o solves SBM-y{—l) and the initial 
condition (uq,Vq) ~ {Uo,Va) do not depend on 7. If furthermore we suppose 

(t/o(A;),yo(fc))e{(0,i),(0,i)}, keZ'\ 

then, for any sequence "fn tending to infinity, we have the convergence in law 

in D{[0,oo), L'p) equipped with the Meyer-Zheng "pseudo-path" topology. Here, U is a standard voter 
process and V ~ \ — U its reciprocal voter process (i.e. opinions 1 and are interchanged) . 

Convention 3.17. In what follows the pair of voter processes constructed in the above theorem will 
be called SBMoo{~l). 



Sketch of Proof. As discussed in the end of Section |1.2[ with the additional assumption on the initial 
conditions, u"^ is a solution to the stepping stone model of Example |1.4| and v'' = 1 — m'^. For 7 
tending to infinity, a well-known result (sec for instance Section 10.3.1 of ^D9l] ) states that the finite 
dimensional distributions of solutions to the stepping stone model converge to those of the standard 
voter process; solutions are bounded and the moments converge as discussed in Section [TTS) Tightness 
in the Meyer-Zheng "pseudo-path" topology follows as for g € (—1, 1). □ 

To understand SBMoo and the voter process in a unified framework let us first summarize. The infinite 
rate symbiotic branching processes SBMoo (g) are the weak limits of SBM^(£<), as 7 — ?► 00, 
• for g € (—1, 1), by Theorem 3.6 



• for g = —1 and Uq -\-Vo = 1, by Theorem |3.16| 

A unified representation can be given with the Poissonian approach developed above if is extended 

to £1 = —1 as 

iy~^{d{vi,V2)) = 5(0,1) (wi, ''^2) ■ 



With the intensities Is{k) defined in (3.11 ) and the Poisson point processes with intensity measure 



A/"' as in (3.10) we can extend Theorem 3.15 as follows 



Theorem 3.18. Suppose g G [—1, 1), (Uq, Vq) G and for g — —1 assume additionally Uq + Vq = 1. 
Then the infinite rate symbiotic branching process {U,V) with initial condition {Uo-,Vo) coincides in 



law with the unique weak solution to [3.13) 
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Note that the additional assumption on the initial condition is not necessary for Equation (3.131 to 
have weak solutions. We believe that also the convergence of SBM^(— 1) to the solutions of (3.13) 
holds without the restriction. 



Proof. For the case g £ (—1, 1) the theorem is nothing else but Theorem 3.15 so that we only need to 
discuss the extension to g = —1. 

Existence of a weak solution to (3.13), for p = —1, can be verified as sketched in the proof of Theorem 



3.13 for g G (—1,1); since the jump measure is finite the proof is simpler since no truncation 
procedure for is needed. 

To identify the weak solutions to (3.13) with SBMoo(— 1) it suffices, by Theorem 3.16 to show that, 
for any weak solution ([/, V) to (3.13), [/ is a voter process and V — I — U . We use two facts: first, 
the jumps preserve the property {Ut{k), Vt(fc)) € {(1, 0), (0, 1)} for alH > 0, fc G Z'' and, secondly, the 
drift and the compensator integral cancel each other. To establish the first, note that th e choi ce of 
implies that always j/i = and y2 — 1 so that the only transitions are (compare with (3.12)) 

-iky 



or, simply. 



Us- 




Vs-{k) 





Us-{k) 
K-(fc) 




The latter follows from the simple computation 



VM-Us{k) 

Us{k)-Vsik) 



(fc) 
(fc) 

AC/,(fc) 



(2/1 



Us- 



l{C/,(/c)=0} + 



(k) 
(k) 



W{{k},d{yi,y2),dr, ds) 



jo AUs{k)l{u,(k)=o} - AVs{k)li^v,{k)=o}ds \ 

{v,(k)=o} ds J 

'jlAUs{k)ds\ 
jlAVs{k)ds 



4v,(fc)=o} 



ds 



for which we used Us{k),Vs{k) G {0, 1} and 

1 



A[/,(fc) + Ay,(fc) = ^(^^■(^■) + ^''(^■)) - ('^^(^) + ^^■(^)) 



lj-fcl=l 



E 



1 

2d 



1 = 0. 



Hence, canceling the compensator integral with the drift shows that Equation (3.13) can be written 
equivalently in the simplified form 



Ut{k) 
Vt{k) 



Uo{k) 
Voik) 



^0 



Us- 



■^f{{k},d{yi,y2),dr, ds). 



Since the configurations only change by a jump and the jumps only switch to 1 and vice versa one 
can already guess that both coordinates are reciprocal voter processes. To make this precise we apply 
Ito's formula to functions of {Ut,Vt) and derive that {U,V) satisfies the martingale problem for the 
standard voter process. It suffices to carry this out for U since we already know that Vt = 1 — Ut , for 
all t > 0. 
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Let US fix a test-function / : {(0, 1)}^ — > M that only depends on finitely many coordinates k E K, 
< oo, and apply Ito's formula to f{Ut) to obtain 



f{Ut) = f{Uo) + J2 



keK 



/(({/,_)«) -/(C/,_) mk},d{y,,y2),dr,ds). 



We denoted again by ?7^'"'-' the configuration that is obtained from the configuration rj fiipping only the 
opinion at site k. Adding and subtracting the compensated integral leads to 

f{Ut)^f{Uo) 



E 



E 

keK 

t 

I. 







/((C/,_)(^^)) -/([/,_) {U-U'){{k},d{yi,y2),dr,ds) 

(fc) [/(([/,_)«) -/(c/,_y 



ds. 



Next, we use that for all s > 0, A: € Z'^ we have Us-{k), Vs-{k) € {0, 1} to obtain 



AU~-'ik) 

V,-{k) 

^ (AVs-{k) 
lAC/,_(fc) 



Us-{k)>0 
V,-{k) > 

: C/,_(fc) = 1 
: Us-{k) = 



= — =1^1 neighbors of the voter at k who have an opinion different than his at time s— } 



= c(fc,C/,). 
Plugging-in, we proved that 



M/ 



fm - fiuo) 



J2c{k,Us) [/((C/,)W) - /(L/, 



keK 



ds 



is a local martingale and since everything is bounded it is, in fact, a martingale. This shows that Ut 



has the generator (1.6 1 of the voter process 



Well-posedness for this martingale problem implies the weak uniqueness statement of the theorem for 
g = -l. □ 

Finally, we want to explain that the extended choice of i^^ is more natural than it appears on first 
view. There are two good reasons. First, going back to Definitions |3.8| and |3. 10| let us see what we get 

for g = —1: 



Van ■ 

e-i-O 



lim - 

e-i-O £ 



:<5(o,i+£) + Yq7^'5(i+e.o) I - S{o^i) + 



- oo5 



(1,0), 



since for completely negatively correlated Brownian motions {B^ , B^) started at (u, v) the exit-measure 



from the first quadrant is 



. Secondly, a more careful look at the density of 



for Q e (—1,1) shows that the mass accumulates at (1,0) and (0,1) since p{q) explodes for g 
tending to —1. More precisely, converges in the vague topology (extended to the completion of M) 
to <5(o.i) -l- cxD(5(]^_g). Unfortunately, both justifications lead to with an additional infinite atom at 
(1, 0). Luckily, the infinite atom at (1, 0) has no impact on the Poissonian equations since the integrand 
of ( 3.13[ ) vanishes if j/2 = and yi ~ 1. We believe that some rigorous work on this observation might 
lead to some interesting results. 

This brief discussion explains the natural unification of the family SBMqo with the voter 
process at its boundary g = —1 and justifies our interpretation of SBMoo(g) as generalized 



voter process, given belov^r Theorem 0.1 
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Appendix A. A Very Rough Primer on Jump SDEs 

Symbiotic branching models are by definition solutions of (possibly infinite) systems of ordinary sto- 
chastic differential equations 

(A.I) dXt = b{Xt)dt + a{Xt)dBt. 

Interestingly, the infinite rate analogues that have been defined so far as solutions to exponential 
martingale problems can be represented as solutions to jump-type stochastic differential equations. 
The most straight- forward generalization of ( A.I[ ) is 



(A.2) dXt = h{Xt)dt + a{Xt)dBt + c{Xt-)dLt 



for a Levy process Lf. The modeling drawback of (A.2 1 is that once Lt has a jump x, then Xt has a 



jump c(A(_)x. If the jumps of the solution process are meant to depend on the jumps of the jump- 
measure in a non-linear way, other concepts are needed. One way to model such processes is to replace 
the jump noise by a general compensated random measure: 

dXt = h{Xt)dt + a{Xt)dBt + c{Xt-,x){Af - N'){dt, dx). 

This notion of jump-type stochastic differential equation is needed for our purposes. Unfortunately, 



the basic jump measure of Theorem O.I has a second order singularity at (I, 0) and a polynomial 
decreasing tail which for q>Q prevents existence of second moments. This causes the general second 
moment integration theory to collapse here and the abstract martingale integration theory with respect 
to compensated random measures comes into play. To guide the reader unfamiliar with those concepts 
we briefly recall some core definitions and concepts. 

First, suppose Af{dt,dx){uj) is a Poisson point measure on [0, oo) x M.'^ with compensator measure A 
on a stochastic basis {Tt)t>o,V), i-c. for all measurable sets A with X{A) < oo, Af{[0,t], A) is a 

Poisson process in t with parameter X{A) such that for disjoint sets Ai, A2 the processes A/'([0, t], Ai) 
are independent. Defining A/"' ([0, t], A) = t\{A), it then follows that the compensated process 

(A.3) iM-N-'m,t],A) :=AA([0,t],A)-AA'([0,t],v4) 

is a martingale. This property motivates the name martingale measure for the random measure Af—Af'. 
Given a predictable integrand H{s, x){uj), defined on the stochastic basis of the driving point process, 
one then aims to define the integral process 

ft p 

H{s,x){U-W){ds,dx) 



of H against the compensated martingale measure Af — Af' via an L^-approximation procedure for 
integrands in the space 

ft 



n = <H -.E 



H{s,xfAf'{ds,dx) 



< 00 



Integrals are then defined as limits of integrals of simple processes against Af — N'. The martingale 
property of the driving measure implies that the stochastic integral itself is a square-integrable mar- 
tingale. For a more detailed introduction we refer the reader for instance to the overview article |B04j 
or [IW8T] . 

Unfortunately, it turns out that for g > our basic jump measure has too heavy tails so that the 
integrands fail to be members of H. Fortunately, abstract martingale theory allows for an integration 
theory with respect to compensated random measures without requiring H ^ %. The rest of this 
section consists of a short summary of the integration theory developed in Section II. Id of |JS03) . To 
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make our life simpler (and this is what we need) we assume that the appearing compensator measure 
is absolutely continuous in t so that the presentation is slightly simplified in contrast to the general 
theory presented in jJS03| . 

Suppose that for a subspace E of M" , A/" is an integer valued random measure (not necessarily Pois- 
sonian) on [0, oo) x E, i.e. a family of measures Af{dt,dx){uj) on ([0, oo) x E^B ® 8) such that 
J\f{{Q} X E){uj) — almost surely, i.e. no jump at time 0, and that Af{-) is an integer. Building 
upon (A. 3 1 the concept of a compensator measure for general random point processes is generalized 



as follows: A/"' is the up to a null set unique (now possibly random) measure such that 
(A.4) * TVf - * A/;' 

is a martingale null at zero for a suitably class of test-functions W. Here, • * -t stands for pathwise 
Lebesgue integration on iJ x [0, t\. As, by assumption, the jump measure Af is integer valued it should 
come as no surprise that Af may be regarded as counting measure for the jumps of an auxiliary 
i^- valued optional process /3t, i.e. 

AA([0,i] X A)(c.)==^l^(A/3,(c.)). 

s<t 

With this notation in hand we can proceed with the abstract definition of the stochastic integral (see 
Definition II. 1.27b) of [JS03 ). Absolute continuity in time of the compensator implies that J\f'{{t} x 
dx){uj) = almost surely so that the quantity W in }JS03| vanishes. The set of possible integrands is 
changed to 



I s<t 



1/2 



< OO 



and the stochastic integral 



Jo J E 



is defined to be the unique (up to indistinguishable) purely discontinuous local martingale Xt such 
that 

(A. 5) AX and iJ(-, A/3.)l{^^ -^q-j are indistinguishable. 

Hence, if N has an atom at (s, x), the stochastic integral H * [N — N') has a jump H{s—, x). Recall 
that by definition a purely discontinuous local martingale is required to be orthogonal to all continuous 
martingales but not to be pathwise everywhere discontinuous. For example, if A^t is a standard Poisson 
process, the compensated process Nt—t is purely discontinuous but far from being pathwise everywhere 
discontinuous. 

The integrability condition for class G is rather unsatisfactory as it involves the jump measure itself 
rather than only its compensator which might be more easy to handle. A characterization of the set 
G is given in Theorem II. 1.33 of |JS03| : it suffices to show that (recall that in our setting W of |JS03| 
vanishes) 

< OO, 



(A.6) 



E 



E 



/ / H'^{s,x)l[\H{s,x)\<i}Af'ids,dx) 
H{s- x)\l{\H{s,x)\>i}Af' {ds, dx) 



< oo. 



showing in particular that G C H. Finally, to motivate the naming "stochastic integral" for the 
abstract local martingale H* {J\f — Af')t, the following property should be mentioned. If the integrand 
is nice, that is, additionally E[|_ff| * A/"/] < oo, then both integrals against Af and the compensator 
measure Af' can be defined pathwise and 

(A.7) H*{Af-Af')t=H*Aft-H*Afi. 
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